Nambu-Jona-Lasinio model description of weakly interacting Bose condensate 
and BEC-BCS crossover in dense QCD-like theories 
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QCD-like theories possess a positively definite fermion determinant at finite baryon chemical potential fi s 
and the lattice simulation can be successfully performed. While the chiral perturbation theories are sufficient to 
describe the Bose condensate at low density, to describe the crossover from Bose-Einstein condensation (BEC) 
to BCS superfluidity at moderate density we should use some fermionic effective model of QCD, such as the 
Nambu-Jona-Lasinio model. In this paper, using two-color two-flavor QCD as an example, we examine how the 
Nambu-Jona-Lasinio model describes the weakly interacting Bose condensate at low density and the BEC-BCS 
crossover at moderate density. Near the quantum phase transition point /jb = m„ (fn* is the mass of pion/diquark 
multiplet), the Ginzburg-Landau free energy at the mean-field level can be reduced to the Gross-Pitaevskii free 
energy describing a weakly repulsive Bose condensate with a diquark-diquark scattering length identical to 
that predicted by the chiral perturbation theories. The Goldstone mode recovers the Bogoliubov excitation in 
weakly interacting Bose condensates. The results of in-medium chiral and diquark condensates predicted by 
chiral perturbation theories are analytically recovered. The BEC-BCS crossover and meson Mott transition at 
moderate baryon chemical potential as well as the beyond-mean-field corrections are studied. Part of our results 
can also be applied to real QCD at finite baryon or isospin chemical potential. 

PACS numbers: ll.10.Wx, 12.38.-t, 25.75.Nq 



I. INTRODUCTION 

A good understanding of quantum chromodynamics (QCD) 
at finite temperature and baryon density is crucial for us to un- 
derstand a wide range of physical phenomena. For instance, 
to understand the evolution of the Universe in the first few 
seconds, one needs the knowledge of QCD phase transition 
at temperature T ~ 200MeV and very small baryon density. 
On the other hand, understanding the physics of neutron stars 
requires the knowledge of QCD at high baryon density and 
very low temperature |l|j. Lattice simulation of QCD at fi- 
nite temperature has been successfully performed in the past 
few decades; however, no successful lattice simulation at high 
baryon density has been done due to the sign problem 0,[3t]: 
The fermion determinant is not positively definite in presence 
of a nonzero baryon chemical potential ji B . 

We thus look for some special cases which have a posi- 
tively definite fermion determinant. One case is QCD at fi- 
nite isospin chemical potential /ui 0, @1, where the ground 
state changes from a pion condensate to a BCS superfluid with 
quark-antiquark condensation with increasing isospin density. 
Another case is the QCD-like theories J6l- tl2ll where quarks 
are in a real or pseudoreal representation of the gauge group, 
including two-color QCD with quarks in the fundamental rep- 
resentation and QCD with quarks in the adjoint representa- 
tion. While these cases do not correspond to the real world, 
they can be simulated on the lattice and may give us some in- 
formation of real QCD at finite baryon density. For all these 
special cases, chiral perturbation theories predict a continuous 
quantum phase transition from the vacuum to the matter phase 
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at baryon or isospin chemical potential equal to the pion mass, 
in contrast to real QCD where the phase transition takes place 
at fi B approximately equal to the nucleon mass. The result- 
ing matter near the quantum phase transition is a dilute Bose 
condensate of diquarks or pions with weakly repulsive interac- 
tions IU3I1 . The equations of state and elementary excitations 
in such matter have been investigated many years ago by Bo- 
goliubov 111411 and Lee, Huang, and Yang 11511 . Bose-Einstein 
condensation (BEC) phenomenon is believed to widely exist 
in dense matter, such as pions and kaons can condense in neu- 
tron star matter if the electron chemical potential exceeds the 
effective mass for pions and kaons IU6l - ll9ll . However, the con- 
densation of pions and kaons in neutron star matter is rather 
complicated due to the meson-nucleon interactions in dense 
nuclear medium. On the other hand, at asymptotically high 
density, perturbative QCD calculations show that the ground 
state is a weakly coupled BCS superfluid with the condensa- 
tion of overlapping Cooper pairs IMII5I I20I - I23II . It is interesting 
that the dense BCS superfluid and the dilute Bose condensate 
have the same symmetry breaking pattern and thus are con- 
tinued with one another. In condensed matter physics, this 
phe nomenon was first discussed by Eagles Il24ll and Leggett 
112511 and is now called BEC-BCS crossover. It has been suc- 
cessfully realized using ultracold fermionic atoms in the past 
few years lual . 

While the lattice simulations of two-color QCD at finite 
baryon chemical potential J27l432ll and QCD at finite isospin 
chemical potential II33I - 136I1 have been successfully performed, 
we still ask for some effective models to link the physics of 
Bose condensate and the BCS superfluidity. The chiral per- 
turbation theories Nt- llll 13711 as well as the linear sigma mod- 
els B38I1 . which describe only the physics of Bose condensate, 
do not meet our purpose. The Nambu-Jona-Lasinio (NJL) 
model [39] with quarks as elementary blocks, which describes 



well the mechanism of chiral symmetry breaking and low en- 
ergy phenomenology of the QCD vacuum, is generally be- 
lieved to work at low and moderate temperatures and densi- 
ties l40m2ll . Recently, this model has been used to describe 
the superfluid transition at finite chemical potentials IU2ll43l - 
DIII53I - D6II for the special cases we are interested in this paper. 
One finds that the critical chemical potential for the superfluid 
transition predicted by the NJL model is indeed equal to the 
pion mass J48l 15311 . and the chiral and diquark condensates 
obtained from the mean-field calculation agree with the re- 
sults from lattice simulations and chiral perturbation theories 
near the quantum phase transition J48l 15311 . The NJL model 
also predicts a BEC-BCS crossover when the chemical po- 
tential increases J50l l54l - D6ll . A natural problem arises: how 
can the fermionic NJL model describe the weakly interact- 
ing Bose condensate near the quantum phase transition? In 
fact, we do not know how the repulsive interactions among 
diquarks or mesons enter in the pure mean-field calculations 
153145511 . In this paper, we will focus on this problem and show 
that the repulsive interaction is indeed properly included even 
in the mean-field calculations. This phenomenon is in fact 
analogous to the BCS description of the molecular conden- 
sation in strongly interacting Fermi gases studied by Leggett 
many years ago 12511 . Fermionic models have been used to de- 
scribe the BEC-BCS crossover in cold Fermi gases by the cold 
atom community. Recently, it has been shown that we can re- 
cover the equation of state of the dilute Bose condensate with 
correct boson-boson scattering length in the strong coupling 
limit, including the Lee-Huang-Yang correction by consider- 
ing the beyond-mean-field corrections [57-59]. In Appendix 
lAl we give a summary of the many -body theoretical approach 
in cold atoms, which is useful for us to understand the theo- 
retical approach and the results of this paper. 

In this paper, using two-color two-flavor QCD as an exam- 
ple and following the theoretical approach of the BEC-BCS 
crossover in cold Fermi gases 0571 15811 . we examine how the 
NJL model describes the weakly interacting Bose conden- 
sate and the BEC-BCS crossover. Near the quantum phase 
transition point /i B - m„, we perform a Ginzburg-Landau 
expansion of the effective potential at the mean-field level, 
and show that the Ginzburg-Landau free energy is essentially 
the Gross-Pitaevskii free energy describing weakly interact- 
ing Bose condensates via a proper redefinition of the conden- 
sate wave function. As a by-product, we obtain a diquark- 
diquark scattering length «dd = m n l(\6nf%) (f„ is the pion 
decay constant) characterizing the repulsive interaction be- 
tween the diquarks, which recovers the tree-level result pre- 
dicted by chiral Lagrangian lla- Ulll . We also show analytically 
that the Goldstone mode takes the same dispersion as the Bo- 
goliubov excitation in weakly interacting Bose condensates, 
which gives a diquark-diquark scattering length identical to 
that in the Gross-Pitaevskii free energy. The mixing between 
the sigma meson and diquarks plays an important role in re- 
covering the Bogoliubov excitation. The results of in-medium 
chiral and diquark condensates predicted by chiral perturba- 
tion theory are analytically recovered. At high density, we 
find the superfluid matter undergoes a BEC-BCS crossover at 
A<b - {ma- 1 m,,) 1 ^ m n =* (1.6 - 2)m„- with m^ being the mass 



of the sigma meson. At /ib - 3m^, we find that the chiral 
symmetry is approximated restored and the spectra of pions 
and sigma meson become nearly degenerated. Well above the 
chemical potential of chiral symmetry restoration, the degen- 
erate pions and sigma meson undergo a Mott transition, where 
they become unstable resonances. Because of the spontaneous 
breaking of baryon number symmetry, mesons can decay into 
quark pairs in the superfluid medium at nonzero momentum. 

The beyond-mean-field corrections are studied. The ther- 
modynamic potential including the Gaussian fluctuations is 
derived. It is shown that the vacuum state I//3 1 < m n is ther- 
modynamically consistent in the Gaussian approximation, i.e., 
all thermodynamic quantities keep vanishing in the regime 
I/^bI < m* even though the beyond-mean-field corrections are 
included. Near the quantum phase transition point, we expand 
the fluctuation contribution to the thermodynamic potential in 
powers of the superfluid order parameter. To leading order, 
the beyond-mean-field correction is quartic and its effect is to 
renormalize the diquark-diquark scattering length. The cor- 
rection to the mean-field result is shown to be proportional to 
m^/f^. Thus, our theoretical approach provides a new way 
to calculate the diquark-diquark or meson-meson scattering 
lengths in the NJL model beyond the mean-field approxima- 
tion. We also find that we can obtain a correct transition tem- 
perature of Bose condensation in the dilute limit, including 
the beyond-mean-field corrections. 

The paper is organized as follows: In Sec. [TT] we derive the 
general effective action of the two-color NJL model at finite 
temperature and density, and determine the model parameters 
via the vacuum phenomenology. In Sec. [TTTJ we investigate the 
properties of dilute Bose condensate near the quantum phase 
transition at the mean-field level. In Sec. [V] the properties 
of matter at high density are discussed. Beyond-mean-field 
corrections are studied in Sec. IIVI We summarize in Sec. [VI] 
Natural units are used throughout. 



II. NJL MODEL OF TWO-COLOR QCD 

Without loss of generality, we study in this paper two-color 
QCD (the number of colors N c - 2) at finite baryon chemical 
potential jUb . For vanishing current quark mass too, two-color 
QCD possesses an enlarged flavor symmetry S\J(2Nf) [Nf is 
the number of flavors], the so-called Pauli-Gursey symmetry 
which connects quarks and antiquarks iq- Tllll . For Nf - 2, 
the flavor symmetry SU(4) is spontaneously broken down to 
Sp(4) driven by a nonzero quark condensate (qq) and there 
arise five Goldstone bosons: three pions and two scalar di- 
quarks. For nonvanishing current quark mass, the flavor sym- 
metry is explicitly broken, resulting in five pseudo-Goldstone 
bosons with a small degenerate mass m„. At finite baryon 
chemical potential /i B , the flavor symmetry SU(2iV/) is ex- 
plicitly broken down to SU l (A ? /)®SU r (A^/)®Ub(1). Further, 
a nonzero diquark condensate {qq) can form at large enough 
chemical potentials and breaks spontaneously the Ub(1) sym- 
metry. In two-color QCD, the scalar diquarks are in fact the 
lightest "baryons," and we expect a baryon superfluid phase 
with (qq) + for |/i B | > m„. 



To construct a NJL model for two-color two-flavor QCD 
with the above flavor symmetry, we consider a contact 
current-current interaction GaYj^^iqy^t^iqy^t^q) where f a 
(a = 1,2,3) are the generators of color SU C (2) and G c is a 



J 



phenomenological coupling constant. After the Fierz trans- 
formation we can obtain an effective NJL Lagrangian density 
with scalar mesons and color singlet scalar diquarks 113311 - 



Xnjl = qiyfdy. - mo)q + G [(qq) 2 + (qiy 5 Tq) 2 + (qiy5T 2 tiqc)(qciy5T2hq)\ , 

I 



(1) 



where q c - Cq T and q c - q T C are the charge conjugate 
spinors with C = ijoji and t; (i = 1,2, 3) are the Pauli matri- 
ces in the flavor space. The four-fermion coupling constants 
for the scalar mesons and diquarks are the same, G - 3G c /4 
15311 . which ensures the enlarged flavor symmetry SU(2iV/) of 
two-color QCD in the chiral limit mo = 0. One can show 
explicitly that there are five Goldstone bosons (three pions 
and two diquarks) driven by a nonzero quark condensate (qq). 
With explicit chiral symmetry broken mo + 0, pions and di- 
quarks are also degenerate, and their mass m n can be deter- 
mined via the standard method for the NJL model |4C 



where we adopt the finite temperature formalism with r = it, 

x - (t, r), and J dx — L dr J d 3 r. The partition function 
can be bosonized after introducing the auxiliary boson fields 



<t(x) - -2Gq(x)q(x), n(x) - -2Gq(x)iy=,Tq(x) (3) 

for mesons and 



4>(x) - -2Gq c (x)iy 5 T 2 t2q(x) 



(4) 



A. Effective action at finite temperature and density 

The partition function of the two-color NJL model (Q]) at 
finite temperature T and baryon chemical potential fis is 



-NJL 



I [ dq] [dq] exp I dx I 



\-Cwl + Y q7oq ) 



(2) 



for diquarks. With the help of the Nambu-Gor'kov represen- 
tation ^ - ( q q c \ the partition function can be written as 

Z N jl = |W] [dW] [da] [dx] [d<ft] [dcf>] exp (-SUa) , (5) 
where the action y[ e ff is given by 



^ = j dx crHx) + n\x) + \4>(x)? _ j ^ J^^- WW) 



with the inverse quark propagator defined as 



G _1 (x,x') = 



y°(-d T + If) + iy • V - Mix) 

-iy54> J (x)T 2 t2 



-iy5<p(x)T 2 t2 



y°(-d T - f ) + iy ■ V - M T (x) 



TVvA \5(x-x). 



(6) 



(7) 



Here Mix) — mo + cr(x) + iy^r ■ n(x). After integrating out the quarks, we can reduce the partition function to Znjl = 
J[dcr] [dn] [d<p'] [d(p] exp - >S e ff [<x, n, <j>\ cf>] J, where the bosonized effective action S e ff is given by 



S e ti[cr, 7t, <p , 4>] - I dx — — 



^!-I T rlnG-'(,,A 
2 



(8) 



Here the trace Tr is taken over color, flavor, spin, Nambu- 
Gor'kov and coordinate (x and x') spaces. The thermody- 
namic potential density of the system is given by Q.(T,jub) - 
- limy^coiT/V) lnZ NJL . 



B. Evaluating the effective action 

The effective action S e g as well as the thermodynamic po- 
tential Q cannot be evaluated exactly in our 3 + 1 dimensional 
case. In this work, we firstly consider the saddle point approx- 
imation, i.e., the mean-field approximation. Then we investi- 
gate the fluctuations around the mean field. 



(I)Mean-field approximation. In this approximation, all 
bosonic auxiliary fields are replaced by their expectation val- 
ues. To this end, we write (cr(x)) = u, (c/>(x)) - A and set 
(7r(x)) - 0. While A can be set to be real, we do not do this 
first in our derivations. We will show in the following that 
all physical results depend only on |A| 2 . The zeroth order or 
mean-field effective action reads 



Q(0) - Y _ 
°eff _ j 



v 2 + |A| 2 i v g-\K) 



(9) 



Here and in the following K - (ia>„, k) with d„ - (2n + \)nT 
being the fermion Matsubara frequency, and Y*k - T 2„ 2k 
with 2k = J<i 3 k/(27r) 3 . The inverse of the Nambu-Gor'kov 
quark propagator §~ x (K) is given by 



(m, + f ) r ° - y k 

-iy 5 A i T 2 t2 



M - iy 5 At 2 t 2 

{iiO n - f )r° - y ■ k 



M 



(10) 



with the effective quark Dirac mass M = mo + v. The mean- 
field thermodynamic potential Qo 
ated as 



,«» 



(T/V)S"' can be evalu- 



'eff 



O.0 



v 2 + |A| 2 
AG 



2N c N f J] [nV(£+) + W(E k )] (11) 



with the definitions of the function "W(E) = E/2 + 
rin(l + e~ E l T ) and the BCS-like quasiparticle dispersions 
£* = V( £ k ±^/ B /2) 2 + |A| 2 where E k = Vk 2 + M 2 . The 
signs + correspond to quasiquark and quasi-antiquark exci- 
tations, respectively. The integral over the quark momentum 
k is divergent at large |k|, and some regularization scheme 
should be adopted. In this paper, we employ a hard three- 
momentum cutoff A. 

The physical values of the variational parameters M (or v) 
and A should be determined by the saddle point condition 



5v 



&Sf>,A] 
8A 



(12) 



-.mi 



which minimizes the mean-field effective action S „. One can 

elf 

show that the saddle point condition is equivalent to the fol- 
lowing Green function relations 

(qq) = £Ti0ii(X), 

K 

(qdysTitiq) = J]Tt[gn(K)iy5T2t2] , (13) 

K 

where the matrix elements of Q are explicitly given by 



0u{K) = 

&22(K) = 

Sn(K) = 

02i{K) = 



(^) 2 -w Akro+ 



■£ 



(iaj,,) 2 - (£ k ) 2 

-iArjfy 
(ico n ) 2 - (£ k ) 2 

-iAh 2 t2 



~^\7o + 



Krs + 



{ico n ) 2 - (E.) 2 



~Z2 A k75 + 



(i^-iEt) 2 ^ ' 

-iAr 2 t2 
(^„) 2 -(£+/' 

-iA T T 2 t 2 



A kT5 



with the help of the massive energy projectors 16] 



Ak 



1 ± 



ro (r ■ k + M) 



E k 



(15) 



Here we have defined the notation £ k = £k + Mb /2 for conve- 
nience. 

(II)Derivative expansion. Next, we consider the fluctuations 
around the mean field, corresponding to the bosonic collective 
excitations. Making the field shifts for the auxiliary fields, 

o-(x) — > u + o-(x), n(_x) — > + n(x), 

cp{x) -> A + (f){x), cj)\x) -» A" 1 ' + (f>\x), (16) 

we can express the total effective action as 



c e(0) ^ f a (o 2 +* 1 + \4>? ^ vo- + A0 + + A" 1 " 

6 e ff = o„+ I dx\ 

eft eff J \ 4G 2G 

- -Trln 1+ I £/xi£(x,jti )£(*!, x') 



(17) 



Here Q(x, x') is the Fourier transformation of 0(ia> n ,k), and 
Y,(x, x') is defined as 



Sfr x ') = ( ~°^ ~ ij5T ' n ^ -'75<P( x ) T 2t2 

\ -iy^(x)j 2 t 2 -o-(x) - iy 5 T T ■ 7C(x) 

x S(x - x). 
With the help of the derivative expansion 

" (-D" +1 



K — L - Tr[^E]", 



(18) 



(19) 



we can calculate the effective action in powers of the fluctua- 
tions cr(x), n(x), 4>{x), <p\x). 

The first order effective action S 1 which includes linear 

eft 

terms of the fluctuations should vanish exactly, since the ex- 
pectation value of the fluctuations should be exactly zero. In 
fact, S 1 can be evaluated as 

eit 



*2 = f dX { 



V 

2G 



1 



-ji(g u +@22) 



o-(x) 



+ -Tr[/r 5 (^iiT + fe-r T )]-^(x) 

cf>Hx) 



A 1 
^ + T^r(iy 5 T 2 t 2 gi 2 ) 



A" 1 " 1 

^ + 2 Tr ('T5T2f2^2l) 



4>(x) 



(20) 



We observe that the coefficients of n(x) is automatically zero 
after taking the trace in Dirac spin space. The coefficients of 
<p{x), <p\x) and <x(x) vanish once the quark propagator takes 
the mean-field form and M, A take the physical values satisfy- 
ing the saddle point condition. Thus, in the present approach, 
the saddle point condition plays a crucial role in having a van- 
ishing linear term in the expansion. 

The quadratic term S^I corresponds to the Gaussian fluctu- 



(ioj n ) 2 - (£+) 2 



-A k y5(14) ations. It reads 



/■ 



8) = . ,/WtrW^Mf 1 
4G 4 



/ 



dx\dxidxiQ{x, x\)T.(x\,X2)G(x2, x 3 )£(x 3 , x') 



(21) 



For the convenience of our investigation in the following, we ing bilinear form 
will use the form of S ^ in the momentum space. After the 
Fourier transformation, it can be written as 



? (2) 

'eff 



-kz. 



\QJQtf + HQ)\ 2 + \<p{Q)\ 2 
2G 



+ ^Tr[£C*DX(-0£CK '+ 02(G)]}- (22) 



where Q — (iv m , q) with v ln - 2nmT being the boson Matsub- 
ara frequency and YjQ = T £ m Zq- Here A(Q) is the Fourier 
transformation of the field A(x), and X(0 is defined as 16111 



£(G) = 



-oiQ) - ijsr ■ ;r(0 



-iy5<f>(Q)T 2 t2 



-iy 5 f(-Q)T2t2 -<t{Q) - iy 5 r l ■ n(Q) 



■ (23) 



(III)Gaussian fluctuations. After taking the trace in Nambu- 
Gor'kov space, we find that «S „! can be written in the follow- 



5 S = ^ZU t(G) * { - Q) crt (e))M(0 

e 
+ ?Z(*i (fi) *2(G> ^(G))N(G) 



( 0(G) 
> + (-G) 

cr(G) 



^2(G) 
^(0 



The matrix M takes the following nondiagonal form 
(ic +U n(Q) n i2(G) 



(24) 



M(G) 



n 13 (G) 
n 23 (G) 
n 31 (G) n 32 (g) ^ + n 33 (0j 



n 2 i(G) 37;+n22(G) 



(25) 

The polarization functions TkjiQ) (i,j = 1,2,3) are one-loop 
susceptibilities composed of the matrix elements the Nambu- 
Gor'kov quark propagator, and can be expressed as 



rin(G) = \Y Tr[g22(K)rgu(P)rh n22(e) = \Yj Tx[gn{K)Tgil{p)V] > 

K K 

niaCfi) = \ Yj Tr [&n{K)VQ n {P)Y] , n 2 i(fi) = \ Y Tr fe(W 21 (P)Il , 
n 33 (G) = \ Yj Tr \8u{K)@u{P) + &22{K)@22(P) + @n(K)g 2l (P) + &2i{K)@ n {P)~\ , 

n 13 (G) = \j] TT[gn(K)rgn(p) + g22(K)rgi2(p)h n3l(e) = \Y Tv[02l{K)gu{P)Y + g2l{K)g2l{p)T] ' 

K K 

n 23 (G) = \ Y Tr ^iiffl r fe( p ) + ^2i(^)r^ 22 ( J P)] , n 32 (0 = I £ Tr Ki,(Z)g 12 (P)r + ^fflfeCT , (26) 



where P = 7T + T = /y5T2?2 and the trace is taken over 
color, flavor and spin spaces. Using the fact that QriiK,^) - 
0u(K,-/Jb) and @ 2 i(K,Hb) = Q\ 2 {K,-^ B ), we can easily 
show that 

n 2 2(0 = n n (-0, n 12 (0 = i4(0, 
n I3 (0 = n^(0 = i4(-0 = n 32 (-G). (27) 

Therefore, only five of the polarization functions are indepen- 
dent. At T = 0, their explicit form is shown in Appendix iBl 
For general case, we can show that II12 °c A 2 and n 13 oc MA. 
Thus, in the normal phase where A = 0, the matrix M recov- 
ers the diagonal form. The off-diagonal elements II13 and II23 



r 



represents the mixing between the sigma meson and diquarks. 
At large chemical potentials where the chiral symmetry is ap- 
proximately restored, M — > mo, this mixing can be safely ne- 
glected. 

On the other hand, the matrix N of the pion sector is diago- 
nal and proportional to the identity matrix, i.e., 



Nij(g) = £j 



^ + n„ ( 



ij = 1,2,3. 



(28) 



This means pions are eigen mesonic excitations even in the 



superfluid phase. The polarization function 11^(0 is given by 

Hr(0 = \Yl T \ gn{K)iy5gn{P)iyi+g2l{K)iy5g22{P)ij5 



- 0n(K)iy 5 g 2 i(P)i75-&2i(K)iy 5 &i2(P)i75 



(29) 



Its explicit form at T — is shown in Appendix iBl We find 
that Yl K (Q) and 1133(0 is different only to a term proportional 
to M 2 . Thus, at high density where (qq) — > 0, the spectra of 
pions and sigma meson become nearly degenerate which rep- 
resents the approximate restoration of chiral symmetry. 
(IV)Goldstone's theorem. The Ub(1) baryon number sym- 
metry is spontaneously broken by the nonzero diquark con- 
densate {qq) in the superfluid phase, resulting in one Gold- 
stone boson. In our model, this is ensured by the condition 
detM(g = 0) = 0. From the explicit form of the polariza- 
tion functions shown in Appendix [B] we find that this condi- 
tion holds if and only if the saddle point condition ( TTZb for v 
and A is satisfied. We thus emphasize that in our theoretical 
framework, the condensates v and A should be determined by 
the saddle point condition, and the beyond-mean-field correc- 
tions are possible only through the thermodynamics, i.e., the 
equations of state. 



C. Vacuum and model parameter fixing 

For a better understanding our derivation in the following, 
it is useful to review the vacuum state at T — /j.b = 0. In the 
vacuum, it is evident that A = and the mean-field effective 
potential Q vac can be evaluated as 



^vac(M) = 



(M - m ) 2 
4G 



■2N c N f Y^E^. (30) 



The physical value of M, denoted by Af«, satisfies the saddle 
point condition dQ. rdc /dM - and minimizes Q. vac . 

The meson and diquark excitations can be obtained from 
tSLp which in the vacuum can be expressed as 



? (2) 

eff 



Q 



o-(-Q)&- 1 (QMQ) 



2 J (2n 

3 

2 



where 



(31) 
respec- 



>i,02 are the real and imaginary parts of 
tively. The inverse propagators in vacuum can be expressed in 
a symmetrical form 14111 

£>r'(e) = ^ + n;(0, 1 = 0-^,4, 

nj(0 = 2iN c N f (Q 2 - ef)I(Q 2 ) 

- 4 ^/0CT (32) 



where e^ = 2M t , e„ — e$ - 0, and the function I(Q 2 ) is 
defined as 



w 



■>-/ 



d 4 K 



1 



(2;r) 4 (Kl - M 2 )(K 2 - M 2 ) ' 



(33) 



with K ± — K + 02. Keeping in mind that M* satisfies the 
saddle point condition, we find that the pions and diquarks are 
Nambu-Goldstone bosons in the chiral limit, corresponding 
to the symmetry breaking pattern SU(4) — >Sp(4). Using the 
gap equation of M„, we find that the masses of mesons and 
diquarks can be determined by the equation 



m {) 



1 



M, 4iGN c N f I(m 2 ) 



+ ef 



(34) 



Since the Q dependence of the function I(Q ) is very weak, 
we find m 2 ~ mo and m 2 ^ =* AM 2 + m 2 . 

Since pions and diquarks are deep bound states, their propa- 
gators can be well approximated by D*„{Q) - -g 2 /(Q 2 -m 2 ) 
with g~y q =* -2iN c NfI(0). The pion decay constant f„ can be 
determined by the matrix element of the axial current, 

= -^Tr f^i [y^sngiK^g^ysT^iK.)] 

= 2N c N fgmM M t Q fl I(Q 2 )S ij . (35) 



1 . Thus, the pion decay constant 
2iN c N f M 2 I(0). (36) 



Here Q(K) = (y"^ - M,) 
can be expressed as 

f 2 

Finally, together with d34T l and d36l l. we recover the well- 
known Gell-Mann-Oakes-Renner relation m 2 f 2 - -mo(qq)o. 
There are three parameters in our model, the current quark 
mass mo, the coupling constant G and the cutoff A. In princi- 
ple they should be determined from the known values of the 
pion mass m K , the pion decay constant f„ and the quark con- 
densate (^)o- Since two-color QCD does not correspond to 
our real world, we get the above values from the empirical 
values /„ =* 93MeV, {uu) * (250MeV) 3 in the N c = 3 case, 
according to the relation f 2 , (qq)o ~ N c . To obtain the model 
parameters, we fix the values of the pion decay constant f„ and 
slightly vary the values of the chiral condensate (qq)o and the 



Set 



A [MeV] GA 2 m [MeV] 



657.9 3.105 

583.6 3.676 

565.8 4.238 

565.4 4.776 



4.90 
5.53 
5.43 
5.11 



(uu)l [MeV] M, [MeV] m n [MeV] 



-217.4 
-209.1 
-210.6 
-215.1 



300 
400 
500 
600 



133.6 
134.0 
134.2 
134.4 



TABLE I: Model parameters (3-momentum cutoff A, coupling con- 
stant G, and current quark mass mo) and related quantities (quark 
condensate (uu)o, constituent quark mass M, and pion mass m n ) for 
the two-flavor two-color NJL model Q. The pion decay constant is 
fixed to be /„ = 75 MeV. 



pion mass m n . Thus, we can obtain different sets of model pa- 
rameters corresponding to different values of effective quark 
mass M, and hence different values of the sigma meson mass 
m lT . Four sets of model parameters are shown in Table. U As 
we will show in the following, the physics near the quantum 
phase transition point pb = m n is not sensitive to different 
model parameter sets, since the low energy dynamics is dom- 
inated by the pseudo-Goldstone bosons (i.e., the diquarks). 
However, at high density, the physics becomes sensitive to dif- 
ferent model parameter sets corresponding to different sigma 
meson masses. The predictions by the chiral perturbation the- 
ories should be recovered in the limit m a /m„ — > co. 



III. DILUTE BOSE CONDENSATE: MEAN FIELD 
THEORY 

Now we begin to study the properties of two-color matter 
at finite baryon density. Without loss of generality, we set 
yt/B > 0. In this section, we study the two-color baryonic mat- 
ter in the dilute limit, which forms near the quantum phase 
transition point p% = m n . Since the diquark condensate is 
vanishingly small near the quantum phase transition point, we 
can make a Ginzburg-Landau expansion for the effective ac- 
tion. As we will see below, this corresponds to the mean-field 
theory of weakly interacting dilute Bose condensates. 



Ginzburg-Landau free energy near the quantum phase 
transition 



Since the diquark condensate A is vanishingly small near 
the quantum phase transition, we can derive the Ginzburg- 
Landau free energy functional VglIA(x)] at T — for the 
order parameter field A(x) = (<p(x)) in the static and long- 
wavelength limit. The general form of Vgl[A(x)] can be writ- 
ten as 



Vgl[A(x)] 



HM 



-5-^2 +K J T ~ ^ 



+a\A{x)\ 2 + l -f5\Mx)f 



(37) 



where the coefficients a,/3,y,6,K should be low energy con- 
stants which depend only on the vacuum properties. The cal- 
culation is somewhat similar to the derivation of Ginzburg- 
Landau free energy of a superconductor from the microscopic 
BCS theory 06211 . but for our case there is a difference in that 
we have another variational parameter, i.e., the effective quark 
mass M which should be a function of |A| 2 determined by the 
saddle point condition. 

(I) The potential terms. In the static and long-wavelength 
limit, the coefficients a,/3 of the potential terms can be ob- 
tained from the effective action S e ff in the mean-field approx- 
imation. At T — 0, the mean-field effective action reads 
*-eff = I dx£lo, where the mean-field thermodynamic poten- 



tial is given by 
Q (|A| 2 ,M) = 



(M - m ) 2 + |A| 2 
AG 



-AyV/^(£+ + £ k ). (38) 



The Ginzburg-Landau coefficients a,/3 can be obtained via a 
Taylor expansion of Qq in terms of |A| 2 , 



Q = Qvac(M s ) + a|A| 2 + ^S|A| 4 



0(|A| 6 ), 



(39) 



where Q vac (M„) is the vacuum contribution which should be 
subtracted. One should keep in mind that the effective quark 
mass M is not a fixed parameter, but a function o/|A| 2 via its 
saddle point condition or gap equation dClo/dM — 0. 

For convenience, we define y = |A| 2 . The Ginzburg-Landau 
coefficient a is defined as 



c/QoCy, M) 



dy 
dQ. {y, M) 



v=() 



dy 

dQiiy, M) 



<9Qo(y, M) dM 



dM dy 



dy 



(40) 



where the indirect derivative term vanishes due to the saddle 
point condition for M. After some simple algebra, we get 

1 
"=4G 



N, 



WlsWZ 



E k 



El 

k k 



^b/ 4 



(41) 



where ££ = -0^ + -^?- We can make the above expression 
more meaningful using the pion mass equation in the same 
three-momentum regularization scheme J4U 14211 , 

We therefore obtain a G-independent result 



(42) 



^NeNfiml-fil)^ 



K 



(£* 2 -m 2 /4)(£* 2 -^ 2 /4)' 



(43) 



From the fact that m n <k 2M, and /? > (see below), we 
see clearly that a second order quantum phase transition takes 
place at exactly p% - m n . Thus, the Ginzburg-Landau free 
energy is meaningful only near the quantum phase transition 
point, i.e., \p% - m n \ <§c m^, and a can be further simplified as 



where the factor J' is defined as 



•7=7*^/2 



K 



(£ k 2 



•m 2 /4) 2 



(44) 



(45) 



The coefficient ft of the quartic term can be evaluated via 
the definition 



P 



d 2 £l Q (y, M) 



dy 2 
d 2 £l (y, M) 



dy 2 



v=0 



v=0 



d 2 Q.()(y, M) dM 



dMdy dy 



(46) 



v=() 
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Notice that the last indirect derivative term does not vanish 
here and will be important for us to obtain a correct diquark- 
diquark scattering length. The derivative dM/dy can be ana- 
lytically derived from the gap equation for M. From the fact 
that dClo/dM - 0, we obtain 



d_ ( dClpiy, M) \ d ( dCl (y,M) \ dM _ 
dy\ dM J + dM\ dM )~dy ~ ' 

Thus, we find 

dM _ d 2 Q. (y, M) l d 2 a Q (y,M) 
~dy ~ dMdy \ dM 2 
Then the practical expression for/? can be written as 

P=/3l+/3 2 , 
where f3\ is the direct derivative term 
3 2 Q (j,M) 



A 



and /?2 is the indirect term 



dy 2 



v=() 



# = - 



d 2 D. (y,M)\ (d 2 Cl (y,M) 



dMdy 



dM 2 



(47) 



(48) 



(49) 



(50) 



(51) 



y=0 



Near the quantum phase transition, all chemical potential 
dependence can be absorbed into the coefficient a, and we 
can set yUe = m n in/3. After a simple algebra, the explicit form 
of f}\ and/?2 can be evaluated as 



A 



4 NcNf ^i^(Et-emJ2)^ 



(52) 



and 



fi 2 = -ll NcNf ^ 

l Z e=± k 



M 



1 



El (E* - enix/2) 2 



2GM t 






k c k / 



(53) 



The G-dependent term mo/(2GM„) in (1531 can be approxi- 
mated as m 2 f 2 /M 2 using the relation m 2 f 2 - -mo(qq). 
(II) The kinetic terms. The kinetic terms in the Ginzburg- 
Landau free energy can be derived from the inverse of the di- 
quark propagator 16211 . In the general case with A + 0, the di- 
quarks are mixed with the sigma meson. However, approach- 
ing the quantum phase transition point, A — » 0, the problem is 
simplified. After the analytical continuation iv,„ — > a> + i0 + , 
the inverse of the diquark propagator in the limit /ie — > rn K can 
be evaluated as 



£> d Vq)=^+n d ( W ,q) 



(54) 



where the polarization function IId(w, q) is given by 
n d (co,q) =N c NfY J 



F* + F* 
c k + c k+q 



(w +yUB ) 2 -(£: 



F* ) 2 



k ■ (k + q) + Mi 



i\ 



e:e*. 

k k+q 



(55) 



In the static and long-wavelength limit (a>, |q| — > 0), the 
coefficients k, 6, y can be determined by the Taylor expansion 
D2 1 (a>,q) - £^'(0,0) - 5a> 2 - koj + yq 2 . Notice that a is 
identical to D d l (0, 0) which is in fact the Thouless criterion 
for the superfiuid transition. On the other hand, keeping in 
mind that D d l {a>, q) can be related to the pion propagator in 
the vacuum, i.e., D d l (u,q) - (l/2)£)* _1 (w + yu B ,q), in the 



static and long-wavelength limit and for /i B 
can well approximate it as J4in 



m„ <s 2M„ we 



£> d V q) - -J [(« + A^b) 2 - q 2 - m 2 ] , 



(56) 



where J' is the same factor defined in d44b . and one can show 
that J' =* g~ 2 l2. We thus find that 5 - y — JT which ensures 
the Lorentz invariance of the vacuum, and k =* 2/jb^T. 



B. From Ginzburg-Landau to Gross-Pitaevskii free energy 

We now show how the Ginzburg-Landau free energy can be 
reduced to the theory describing weakly repulsive Bose con- 
densates, i.e., the Gross-Pitaevskii free energy 11631 16411 . 

(I) Nonrelativistic version. First, since the Bose condensed 
matter is indeed dilute, let us consider the nonrelativistic ver- 
sion, where a> «: m„ and the kinetic term oc d 2 /dr 2 is ne- 
glected. To this end, we define the nonrelativistic chemical 
potential H& for diquarks, pd = //b - m n , and further simplify 
the coefficient a as 



a - -fi d {2m„J). 



(57) 



Then the Ginzburg-Landau free energy can be reduced to the 
Gross-Pitaevskii free energy of a dilute repulsive Bose gas, if 
we define a new condensate wave function y(jc) as 



\|/(x) = ^2m lr J'A(x). 
The resulting Gross-Pitaevskii free energy is given by 



(58) 



VopfvW] 



/dx v' (x) 1 1 
L \dr 2m n ) 

-Pd\v(x)\ 2 + -go\y(x)\ 4 



(59) 



where go = ^Odd/nt K . The repulsive diquark-diquark interac- 
tion is characterized by a positive scattering length an defined 



add 



\(mm n 



J 



-2 



(60) 



Set 


1 


2 


3 


4 


Odd according to d60t Dm" 1 ] 


0.0631 


0.0635 


0.0637 


0.0639 


add according to d63ll [m~ l ] 


0.0624 


0.0628 


0.0630 


0.0633 



TABLE II: The values of diquark-diquark scattering length add (in 
units of m^ 1 ) for different model parameter sets. 
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Keep in mind that the scattering length obtained here is at the 
mean-field level. We will discuss the possible beyond-mean- 
field corrections in Sec. [V] Thus, for a dilute medium with 
density n satisfying nah, <sc 1, the system is indeed a weakly 
interacting Bose condensate J13l4l5ll . 

(II) Diquark-diquark scattering length. Even though we 
have shown that the Ginzburg-Landau free energy is indeed a 
Gross-Pitaevskii version near the quantum phase transition, a 
key problem is whether the obtained diquark-diquark scatter- 
ing length Add is quantitatively correct. A numerical calcula- 
tion for ( l60t is straightforward. The obtained values of «dd for 
the four model parameter sets are shown in Table [TTJ We can 
also give an analytical expression based on the formula of the 
pion decay constant in the three-momentum cutoff scheme, 



k £ k 



(61) 



According to the fact that m n <K 2M„ B and Sf can be well 
approximated as 

a a ft _ (2/, 2 /M,) 2 __ f n m\ 



J 



Ml m\fllMl+\fl 4M$' 

f 2 

2M 1 ' 



(62) 



Thus, the diquark-diquark scattering length a^ in the limit 
m jr /(2M t ) — > is related only to the pion mass and decay 
constant, 



«dd 



16^/, 2 



(63) 



The values of «dd for the four model parameter sets accord- 
ing to the above expression are also listed in Table [TT] The 
errors are always about 1 % comparing with the exact numer- 
ical results, which means that the expression (|63l is a good 
approximation for the diquark-diquark scattering length. The 
error should come from the finite value of m K /(2M*). We can 
obtain a correction in powers ofm n /(2M*) J65I1 . but it is obvi- 
ously small, and its explicit form is not shown here. 

The result «dd x m n is universal for the scattering lengths 
of the pseudo-Goldstone bosons. Eventhough the SU(4) flavor 
symmetry is explicitly broken in presence of a nonzero quark 
mass, a descrete symmetry (f>\,(f>2 *-* ^1,^2 holds exactly for 
arbitrary quark mass. This also means that the partition func- 
tion of two-color QCD has a descrete symmetry /i^ <-> /j.i J^]. 
Because of this descrete symmetry of two-color QCD, the an- 
alytical expression d63l of «dd (which is in fact the diquark- 
diquark scattering length in the B - 2 channel) should be iden- 
tical to the pion-pion scattering length at tree level in the 1-2 
channel which was first obtained by Weinberg many years ago 
116611 . Therefore, the mean-field theory can describe not only 
the quantum phase transition to a dilute diquark condensate 
but also the effect of repulsive diquark-diquark interaction. 

(Ill) Equations of state. The mean-field equations of state 
of the dilute diquark condensate are thus determined by the 
Gross-Pitaevskii free energy ( 1591 . Minimizing Vgp W-*)] with 



respect to a uniform condensate \|/, we find the physical mini- 
mum is given by 



IVof = 



go' 



(64) 



and the baryon density is n - |\|/ol . Using the thermodynamic 
relations, we therefore get the well-known results for the pres- 
sure P, the energy density S and the chemical potential //b in 
terms of the baryon density n, 

27ra d d 2 

P(n) - n : , 

m„ 

cv \ 27radd 2 

o(n) = m„n H n , 

m„ 

, . 47Tfl d d 

yU B (n) = m K + n, 



(65) 



which were first obtained by Bogoliubov many years ago \U 
We can examine the above results through a direct numerical 
calculation with the mean-field thermodynamic potential. The 
pressure is given by P — -(Qo - Q V ac) and the baryon density 
reads n - -d£lo/dfiB- In FigQ]we show the numerical results 
for the pressure and the chemical potential as functions of the 
density for the four model parameter sets. At low enough 
density, the equations of state are indeed consistent with the 
results J65] l with the scattering length given by d60l ). It is evi- 
dent that the results at low density are not sensitive to different 
model parameter sets, since the physics at low density should 
be dominated by the pseudo-Goldstone bosons. 

In fact, we can derive the equations of state d65l l analyti- 
cally from the mean-field thermodynamic potential Oq. For 
example, the baryon number density reads 






[E k (E k +^) £ k + (£ k + +4 + )J 



(66) 



Near the quantum phase transition point and to leading order 
of |A| 2 , we obtain 



n ^ V ; j; 



|A| 2 



|A| 2 



(Ef - mJ2) 2 (El + mJ2) 2 



2m„J\A\ 2 = | Vo 



(67) 



only 

the Lee-Huang-Yang corrections IU5I1 which are proportional 
to (na^) 1 ^ 2 are absent in the equations of state. As we have 
shown in Appendix lAl to obtain such corrections, it is nec- 
essary to go beyond the mean field, and a beyond-mean-field 
correction to the scattering length «dd is also possible ll57ll58ll . 
(IV) Relativistic version. We can also consider a relativistic 
version of the Gross-Pitaevskii free energy via defining the 
condensate wave function 



<D(x) = ^A(x). 



(68) 



In this case, the Ginzburg-Landau free energy is reduced to a 
relativistic version of the Gross-Pitaevskii free energy, 



V RG p[<I>( 



x)] - I dx <£> f (x) I 



8t 2 



^'Wh— + 2/"b— -v z |ow 



'5t 



+( m 2_ / ,2 ) | OW |2 + _| <DW |4 



(69) 



The self-interacting coupling A — fij' 2 is now dimensionless 
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and can be approximated by A a m*/f%. For realistic values 
of m^ and /^, we find A ~ 0(1). In this sense, the Bose con- 
densate is not weakly interacting, except for the low density 
limit na 3 dd «c 1 . One should keep in mind that this result can- 
not be applied to high density, since it is valid only near the 
quantum phase transition point. 



J 



0.05 




n / (f m ) 




n / (f m ) 



FIG. 1 : The baryon chemical potential //b and the pressure P as functions of the baryon density n for different model parameter sets. The solid 
lines correspond to the direct mean-field calculation and the dashed lines are given by j65l l. 
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C. Bogoliubov excitation in a dilute diquark condensate 

An ideal Bose-Einstein condensate is not a superfluid. In 
presence of weakly repulsive interactions among the bosons, 
a Goldstone mode which has a linear dispersion in the low 
energy limit appears, and the condensate becomes a super- 
fluid according to Landau's criterion minq[w(q)/|q|] > 0. The 
Goldstone mode which is also called the Bog oliubov mode 
here should have a dispersion given by IU3l4l5ll 



I_q2_/_q2_ 8;ra dd w \ 
V 2m„ \2m K m^ j 



|q| « m n . 



(70) 



Since the Gross-Pitaevskii free energy obtained above is at 
the classical level, to study the bosonic collective excitations 
we should consider the fluctuations around the mean field 
0671 - 16911 . The propagator of the bosonic collective modes is 
given by M~ l (Q) and N~'(0- The Bogoliubov mode cor- 
responds to the lowest excitation obtained from the equa- 
tion detM(<y, q) = 0. With the explicit form of the ma- 
trix elements of M in the superfluid phase, we can analyt- 
ically show that detM(0, 0) = which ensures the Gold- 
stone's theorem. In fact, for (u>, q) = (0,0), we find that 
detM = (Mfj - |M 12 | 2 )M 33 +2|M 13 | 2 (|Mi 2 | -M„). Using the 
saddle point condition for A, we can show that Mn(0, 0) = 
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|Mi2(0, 0)| and hence the Goldstone's theorem holds in the su- 
perfluid phase. Further, we may obtain an analytical expres- 
sion of the velocity of the Bogoliubov mode via a Taylor ex- 
pansion for M{cj, q) around (a>, q) = (0, 0) like those done in 
067I469I1 . Such a calculation for our case is more complicated 
due to the mixing between the sigma meson and diquarks, and 
it cannot give the full dispersion (|70l i. 

On the other hand, since A — > near the quantum phase 
transition point, we can expand the matrix elements of M in 
powers of |A| 2 . The advantage of such an expansion is that 
it cannot only give the full dispersion (TTOb but also link the 
meson properties in the vacuum. Formally, we can write down 
the following expansions: 

Mn(<u, q) = D d l (co, q) + \A\ 2 A(co, q) + 0([A| 4 ), 



M 22 (w,q) 

M 12 (w,q) = M 21 (w,q) 



D- A \-u, q) + \A\ 2 A(-oj, q) + 0(|A| 4 ), 



A 2 B(a>, q) 



■ 0(|A| 4 ), 



ML(a>,q) = AH(co,q) + 0(|A| J ), 



M 13 (w,q) 

M 23 (w,q) 

M 33 (w,q) = £);- 1 («,q) + 0(|A| 2 ). 



M; 2 (w,q) = A'H(-u),q) + <9(|A| J ), 



(71) 



Notice that the effective quark mass M is regarded as a func- 
tion of |A| 2 as we have done in deriving the Ginzburg-Landau 
free energy. Since we are interested in the dispersion in the 
low energy limit, i.e., w, |q| <K m„, we can approximate 
the coefficients of the leading order terms as their values at 
(w,q) = (0,0). That is, 



A(w, q) =* A(-cj, q) =* A(0, 0) = A , 
B(o),q)^B(0,0) = B , 

H(u, q) =* H(-co, q) =* H(0, 0) s H . 



(72) 



Further, since m a » m K , we can approximate the inverse 
sigma propagator D* a T l (a>, q) as its value at (a>, q) = (0,0). 
Therefore, the dispersion of the Goldstone mode in the low 
energy limit can be determined by the following equation: 



det 



(£> d Vq) + |A| 2 A 
A t2 B 
A'7/o 



A 2 B AH 

D d I (-w,q) + |A| 2 A A^H 

AH Dr'(0,0) 



(73) 



r 



Now we can link the coefficients Ao, Bo, ^o and D a l {Q, 0) 
to the derivatives of the mean-field thermodynamic potential 
Q.q and its Ginzburg-Landau coefficients. Firstly, using the 
explicit form of Mi 2 , we find that 



|M 12 (0,0)| = |A| 2 /?i 
Secondly, using the fact that 



B Q =B l 



M„(0,0)-|M 12 (0,0)| = 



dQ 

dJA] 2 



(74) 



(75) 



and together with the definition for A(w, q), 
dMnCy.Af) 



A(w,q) = 



dy 
dM u (y,M) 



v=() 



dy 



dM u (y,M)dM 



8M 



dy 



(76) 



we find the following exact relation: 



A =B + B =B + Bi. (77) 

On the other hand, we have the following relations for Hq and 



or 1 (o,o), 



H n = 



©r (0,0) = 



d 2 Q (y, M) 



dMdy 
d 2 n (y, M) 



dM 2 



(78) 



One can check the above results from the explicit forms of 
Mi 3 and M33 in Appendix|B]directly. Thus we have 



«s 



O^-HO.O) 



= ys 2 . 



(79) 



According to the above relations, Eq. (l73l l can be reduced 



to 



3yS 2 |A| 4 + 2B\A\ 2 [D d : (w, q) + D^-CO, q)] 
+ O d 1 (w,q)£) d 1 (- W ,q) = 0. 



(80) 



It is evident that only the coefficient B appears in the final 
equation. Further, in the nonrelativistic limit a>, |q| «; m„ and 
near the quantum phase transition point, D d ' (a>, q) can be ap- 
proximated as 



£> d '(w,q) * -2m n j[co-^- + yu d 



(81) 
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Together with the mean-field results for the chemical potential 
j^d - golVol 2 = y6|A| 2 /(2m ff JT) and for the baryon density n — 
|\|/o| 2 , we finally get the Bogoliubov dispersion (l70l i. 

We should emphasize that the mixing between the sigma 
meson and the diquarks, denoted by the terms AHq and A ' Hq, 



plays an important role in recovering the correct Bogoliubov 
dispersion. Even though we do get this dispersion, we find the 
procedure is quite different to the standard theory of weakly 
interacting Bose gas IU3l U4l l63l 16411 ■ There, the elementary 
excitation is given only by the diquark-diquark sectors, i.e., 



J 



det 



Mn(G) M 12 (<2) 
M 2 i(0 M 22 (0 



det 



" w + £r _ W + 2gol Vol 2 



1m n 



golYol' 



golVol 



(82) 



But in our case, we cannot get the correct Bogoliubov excita- 
tion if we simply set Hq - and consider only the diquark- 
diquark sector. In fact, this requires Aq - 2Bq - 2B which is 
not true in our case. 

One can also check how the momentum dependence of 
A,B,H and JTj modifies the dispersion. This needs direct 
numerical solution of the equation det M(w, q) = 0. We have 
examined that for \pg - m K \ up to 0.0 lra^, the numerical re- 
sult agrees well with the Bogoliubov formula (T70b . However, 
at higher density, a significant deviation is observed. This is 
in fact a signature of BEC-BCS crossover which will be dis- 
cussed in Sec.lIVI 



D. In-medium chiral condensate 

Up to now we have studied the properties of the dilute Bose 
condensate induced by a small diquark condensate {qq). The 
chiral condensate {qq) will be modified in the medium. In 
such a dilute Bose condensate, we can study the response of 
the chiral condensate to the baryon density n. 

To this end, we expand the effective quark mass M in terms 
of y - |A| 2 . We have 



M-M t = 



dM 
~dy~ 



v=0 



y + 0(y 2 ) 



(83) 



The expansion coefficient can be approximated as 



dM 
~dy~ 



VllM, 



v=0 



mlfJMi+Afl 



1 



2M t 



1+0 



4-Mz 



(84) 



Using the definition of the effective quark mass, M — mo - 
2G(qq), we find that 



(qq)„ 

(qqh 



1 



|A| 2 



4G(qq) M» 



1 



|A| 2 
2M 2 ' 



(85) 



Since the baryon number density reads n = |yo| 2 = 2m JT J'\A\ 2 , 
using the fact that J' =* f 2 /(2M 2 ), we obtain to leading order 



(qq),, 
(qqh 



1 



Vlm n ' 



(86) 



This formula is in fact a two-color analogue of the density 
dependence of the chiral condensate in the N c - 3 case, where 
wehaveSIH 



(qq)n 
(qqh 



1 



fn m n 



(87) 



with 2^n being the pion-nucleon sigma term. In Figf2j we 
show the numerical results via solving the mean-field gap 
equations. One finds that the chiral condensate has a perfect 
linear behavior at low density. For large value of M» ( and 
hence the sigma meson mass m CT ), the linear behavior persists 
even at higher density. 

In fact, the Eq. d86b can be obtained in a model indepen- 
dent way. Applying the Hellmann-Feynman theorem to a di- 
lute diquark gas with energy density S(n) given by J65l ), we 
can obtain d86*l l directly. According to the Hellmann-Feynman 
theorem, we have 



2m ((qq)„ - (qqh) = m 



dS 

d?no 



(88) 



The derivative d&/dmo can be evaluated via the chain rule 
dS/dmo - (dS/dntnXdm^/dmo). Together with the Gell- 
Mann-Oakes-Renner relation m 2 f 2 - -mo(qq)o and the fact 
that da^/drnj, =* aalm^, we can obtain to leading order Eq. 
( 1861 ). Beyond the leading order, we find the correction of or- 
der 0(n 2 ) vanishes. Thus, the next-to-leading order correction 
should be 0(n 5 ^ 2 ) coming from the Lee-Huang-Yang correc- 
tion to the equation of state i72ll . 

Finally, we can show analytically that the "chiral rotation" 
behavior J4T- [llll predicted by the chiral perturbation theories 
is valid in the NJL model near the quantum phase transition. 
In the chiral perturbation theories, the chemical potential de- 
pendence of the chiral and diquark condensates can be analyt- 
ically expressed as 



(qq), 



en 



(qqh 



/4 



(qq)m 
(qqh 



/4 



(89) 



Near the phase transition point, we can expand the above for- 
mula in powers of p^ — ps - m n . To leading order, we have 



(qq)m 
(qqh 



~ 1 



2pd 




(90) 
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Using the mean-field result ( f64-b for the chemical potential fi^, 
one can easily check that the above relations are also valid in 
our NJL model. 



E. Chiral limit 

In the above studies we focused on the "physical point" 
where mo + 0. In the final part of this section, we briefly 
discuss the chiral limit with mo = 0. 

We may naively expect that the results at mo ^0 can be 
directly generalized to the chiral limit via setting m n - 0. The 
ground state is a noninteracting Bose condensate of massless 
diquarks, since m n - and add = 0. However, this cannot be 
true since many divergences develop due to the vanishing pion 
mass. In fact, the conclusion of second order phase transition 
is not correct since the Ginzburg-Landau coefficient /3 van- 
ishes. Instead, the superfluid phase transition is of strongly 



first order in the chiral limit (4 

In the chiral limit, the effective action in the vacuum should 
depend only on the combination a 2 +n 2 + \<p\ 2 due to the ex- 
act flavor symmetry SU(4) =* SO(6). The vacuum is chosen 
to be associated with a nonzero chiral condensate (cr) without 
loss of generality. At zero and at finite chemical potential, the 
thermodynamic potential Qq(M, |A|) has two minima locating 
at (M, | A|) = (a, 0) and (M, |A|) = (0, b). At zero chemical po- 
tential, these two minima are degenerate due to the exact fla- 
vor symmetry. However, at nonzero chemical potential (even 
arbitrarily small), the minimum (0, b) has the lowest free en- 
ergy. Analytically, we can show that b — > M„ at ^ — + . 
This means the superfluid phase transition in the chiral limit 
is of strongly first order, and takes place at arbitrarily small 
chemical potential. Since the effective quark mass M keeps 
vanishing in the superfluid phase, a low density Bose conden- 
sate does not exist in the chiral limit. 
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FIG. 2: The ratio R„ = (qq)„/(qq)o as a function of n/(f^m n ) for different model parameter sets. The dashed line is the linear behavior given 
by®. 
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IV. MATTER AT HIGH DENSITY: BEC-BCS CROSSOVER 
AND MOTT TRANSITION 

The investigations in Sec. [Til] are restricted near the quan- 
tum phase transition point // B = m K . Generally the state of 



matter at high density should not be a relativistic Bose con- 
densate described by J69] i. In fact, perturbative QCD calcu- 
lations show that the matter is a weakly coupled BCS super- 



fluid at asymptotic density 1201 - 12311 . In this section, we will 
discuss the evolution of the superfluid matter as the baryon 
density increases from the NJL model point of view. While 



J 
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some results presented in the following have been published 
elsewhere Id3TJ56I1 . we will still show them for the sake of 
completeness. 





FIG. 3 : The chiral and diquark condensates (in units of (qq)o) as functions of the baryon chemical potential(in units of m n ) for different model 
parameter sets. 
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A. Chiral and diquark condensates 

The numerical results for the chiral condensate (qq) and di- 
quark condensate (qq) are shown in Fig[3] As a comparison, 
we also show the analytical result ( T89b predicted by chiral per- 
turbation theories. While the behavior of the chiral condensate 
is in good agreement with the chiral perturbation theories, the 
diquark condensate deviates significantly from the result (T89b 
for small values of M„ . 

This deviation can be understood from the fact that the chi- 
ral perturbation theories correspond to the nonlinear sigma 
model limit m a — » oo. For finite value of m a , one should 
consider the 0(6) linear sigma model 15611 

Xlsm = ^tp) 2 - imV + \ V - Ha- (91) 

where <p - (cr, tt, 0!,0 2 ) and in 2 < 0. The model parameters 
m 2 ,A,H can be determined from the vacuum phenomenology. 



In this model, we can show that the chiral and diquark con- 
densates are given by 148115611 

(qq)^ m 2 n 



(qqh 

(qq)o 



I ml fil - mL 
l--i+2^ 1 



(92) 



In the nonlinear sigma model limit m^ — > oo, the above re- 
sults are indeed reduced to the result ( |89l predicted by chiral 
perturbation theories. However, for finite values of »v, the re- 
sults can be significantly different from d89b at large chemical 
potential. 



B. BEC-BCS crossover 

While the Ginzburg-Landau free energy can be reduced to 
the Gross-Pitaevskii free energy near the quantum phase tran- 
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sition point, it is not the case at arbitrary /^b- 

When jU B increases, we find that the fermionic excitation 
spectra E^ undergo a characteristic change. Near the quan- 
tum phase transition fis = m* they are nearly degenerate since 
m n <K 2M t and their minima are located at |k| = 0. However, 
at very large yUB the minimum of E^ moves to |k| ^ jUb/2 since 
M — > niQ. Meanwhile the excitation energy of the antifermion 
excitations become much larger than that of the fermion ex- 
citations and can be neglected. This characteristic change of 
the fermionic excitation spectra takes place when the mini- 
mum of the lowest band excitation E^ moves from |k| = 
to |k| * OJM HI MM I22H21, i.e., Mb/2 = M(// B )& 
A schematic plot of this characteristic change is shown in 
FigH The equation ^b/2 = M(fiB) defines the so-called 
crossover point jUb = Mo which can be numerically determined 
by the mean-field gap equations. The numerical results of the 
crossover chemical potential //o for the four model parameter 
sets are shown in TableHITl For reasonable parameter sets, the 
crossover chemical potential is in the range (1.6 - T)m n . 





FIG. 4: A schematic plot of the fermionic excitation spectrum in the 
BEC state (left) and the BCS state (right). 



In fact, an analytical expression for //o can be achieved 
according to the fact that the chiral rotation behavior 
(w)a<b/^W)o - m ^//^B * s St ^ valid in the NJL model at large 
chemical potentials as shown in Fig[3] We obtain [54] 



Ho 

2 



2 

-fM. 

Mo 



Ho * (2MX) 



2x1/3 



(93) 



Using the fact that mo- =* 2Af«, we find that ho can be expressed 

as 



Mo_ 



1/3 



(94) 



Thus, in the nonlinear sigma model limit m cr jm n — > oo, there 
should be no BEC-BCS crossover. On the other hand, this 
means the physical prediction power of the chiral perturbation 
theories is restricted near the quantum phase transition point. 
The fermionic excitation gap A ex (as shown in Figj4]i, de- 
fined as the minimum of the fermionic excitation energy, i.e., 
A ex = mini i {£' k , E^}, can be evaluated as 



y/iM^W 



|Ap 
|A| 



Hb < Ho 
Mb > Mo- 



(95) 



It is evident that the fermionic excitation gap is equal to the su- 
perfluid order parameter only in the BCS regime. This is sim- 
ilar to the BEC-BCS crossover in nonrelativistic systems J67I1 . 
and we find that the corresponding fermion chemical potential 
fi can be defined as fi - Hb/2 - M. The numerical results of 
the fermionic excitation gap A ex for different model parameter 
sets are shown in Fig [5] We find that for a wide range of the 
baryon chemical potential, it is of order <9(M,). The fermionic 
excitation gap is equal to the pairing gap |A| only at the BCS 
side of the crossover, and exhibits a minimum at the quantum 
phase transition point. 

On the other hand, the momentum distributions of quarks 
(denoted by «(k)) and antiquarks (denoted by h(k)) can be 
evaluated using the quark Green function Q\ \ (K). We obtain 



■CB-ifl-fl. 



for quarks, 
for antiquarks. 



(96) 



The numerical results for «(k) and «(k) (for model parameter 
set 1) are shown in Fig(6] Near the quantum phase transition 
point, the quark momentum distribution «(k) is a very smooth 
function in the whole momentum space. In the opposite limit, 
i.e., at large chemical potentials, it approaches unity at |k| = 
and decreases rapidly around the effective "Fermi surface" at 
|k| =* \n\. For the antiquarks, we find that the momentum dis- 
tribution «(k) exhibits a nonmonotonous behavior: it is sup- 
pressed at both low and high densities and is visible only at 
moderate chemical potentials. However, even at very large 
chemical potentials, e.g., yU B = 10m„, the momentum distribu- 
tion n(k) does not approach the standard BCS behavior, which 
means the dense matter is not a weakly coupled BCS super- 
fluid for a wide range of the baryon chemical potential. In 
Fig|7] we show the ratio |A|//z up to ^ =* \0m n . It is clear 
that the ratio is not small even at large chemical potentials. At 
Mb - IO/Mjj, it is about 0.5, which means the dense matter is 
still a strongly coupled BCS superfluid. 

The Goldstone mode also undergoes a characteristic change 
in the BEC-BCS crossover. Near the quantum phase transition 
point, i.e., in the dilute limit, the Goldstone mode recovers 
the Bogoliubov excitation of weakly interacting Bose conden- 
sates. In the opposite limit, we expect the Goldstone mode ap- 
proaches the Anderson-Bogoliubov mode of a weakly coupled 
BCS superfluid, which takes a dispersion w(q) = |q|/ v3 up to 
the two-particle continuum to =* 2|A|. In fact, at large chem- 
ical potentials, we can safely neglect the mixing between the 
sigma meson and diquarks. The Goldstone boson dispersion 



Set 


1 


2 


3 


4 


Crossover chemical potential fio [raj 


1.65 


1.81 


1.95 


2.07 



TABLE III: The crossover chemical potential /j (in units of m n ) for 
different model parameter sets. 
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is thus determined by the equation 



jMu(fi) M 12 (g)\_ 
M 2I (0 M 22 (0 I a 



(97) 



The problem is totally the same as that has been investigated 
in J68ll69ll . Therefore, at very large chemical potentials where 
| A|//y becomes small enough, the Goldstone mode recovers the 
Anderson-Bogoliubov mode of a weakly coupled BCS super- 



fluid. 

Finally, we should emphasize that the existence of a smooth 
crossover from the Bose condensate to the BCS superfluid de- 
pends on whether there exists a deconfinement phase transi- 
tion at finite jUb IBU l32l 18111 and where it takes place. Recent 
lattice calculation predicts a deconfinement crossover which 
occurs at a baryon chemical potential larger than that of the 
BEC-BCS crossover fH. 
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FIG. 5: The fermionic excitation gap A cx (in units of M,) as a function of the baryon chemical potential (in units of m„) for different model 
parameter sets. The effective quark mass M and the pairing gap |A| are also shown by dashed and dash-dotted lines, respectively. 
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C. Chiral restoration and meson Mott transition 



As in real QCD with two quark flavors, we expect the 
chiral symmetry is restored and the spectra of sigma me- 
son and pions become degenerate at high density I82I1 . 
For the two-flavor case and with vanishing mo, the residue 
SU l (2)®SU r (2)®Ub(1) symmetry group at // B + is spon- 
taneously broken down to SpL(2)®SpR(2)in the superfluid 
medium with nonzero {qq) and (qq), resulting in one Gold- 



stone boson. For small nonzero mo, we expect the spectra 
of sigma meson and pions become approximately degenerate 
when the in-medium chiral condensate {qq) becomes small 
enough. 

In fact, according to the result {qq)„/{qq)o - 1 - nj(2f^m 7I ) 
at low density, we can roughly expect that the chiral symme- 
try is approximately restored at n ~ 2f%m n . From the chemi- 
cal potential dependence of the chiral condensate {qq) shown 
in Fig J3j we find that it becomes smaller and smaller as the 
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FIG. 6: The momentum distributions for quarks (upper panel) and antiquarks (lower panel) for various values of /ub- The momentum is scaled 
by|//|=|/i B /2-M|. 
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FIG. 7: The ratio of the pairing gap |A| to the effective fermionic chemical potential /i = /jb/2 - M as a function of the baryon chemical 
potential (in units of m n ) for different model parameter sets. The divergent point corresponds io//b = y"o, i-e., the BEC-BCS crossover point. 



density increases. As a result, we should have nearly degen- 
erate spectra for the sigma meson and pions. To show this we 
need the explicit form of the matrix M(Q) and N(0 given in 
Appendix iBl Since Mb,M32 <* MA, at high density where 
(qq) — > 0, they can be safely neglected and the sigma meson 
decouples from the diquarks. The propagator of the sigma 
meson is then given by M^ 1 (Q). From the explicit form of the 
polarization functions Yla-(Q) = ^33(Q) and Tl^iQ), we can 



see that the inverse propagators of the sigma meson and pi- 
ons differ from each other in a term proportional to M 2 . Thus 
at high density their spectra are nearly degenerate, and their 
masses are given by the equation 



l-2Gn ff (w,0) = 0. 



(98) 



Using the mean-field gap equation for A, we find the solution 
is <d — /is, which means the meson masses are equal to /j.b 
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at large chemical potentials. In Fig|8] we show the chemical 
potential dependence of the meson and diquark mass spec- 
tra determined at zero momentum. We find from the meson 
spectra that the chiral symmetry is approximately restored at 
/iB — 3m,r, corresponding to n =* 3.5f^m K . It is interesting 
that near the quantum phase transition point /iB = m n the mix- 
ing between the sigma meson and diquarks is very strong and 
makes the sigma meson lost its way. Since it is continuous 
with the antidiquark mode in the normal phase, it is also called 
the "antidiquark" mode in the superfluid phase 11531 18311 . The 
"sigma meson," which is continuous with the sigma meson in 
the normal phase, is in fact the Higgs mode of the BCS super- 
fluid with a mass 2A at high density. 



FIG. 8: The mass spectra of mesons and diquarks (in units of m n ) as 
functions of the baryon chemical potential (in units of m n ) for model 
parameter set 1. For other model parameter sets, the mass of the 
heaviest mode is changed but others are almost the same. 
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FIG. 9: The two-particle continua ojq q and a> qq (in units of M,) as functions of the baryon chemical potential (in units of m„) for different model 
parameter sets. The degenerate mass of pions and sigma meson is shown by dashed line. 
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Even though the deconfinement transition or crossover meson Mott transition associated with the chiral restoration 
which corresponds to the gauge field sector cannot be de- I84486I1 . The meson Mott transition is defined as the point 
scribed in the NJL model, we can on the other hand study the where the meson energy becomes larger than the two-particle 
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continuum a) qq for the decay process n — > qq at zero mo- 
mentum, which means the mesons are no longer bound states. 
The two-particle continuum a> qq is different at the BEC and 
the BCS sides. From the explicit form of Yl„(Q), we find that 



<Jqq 



|A|+ ^(M+^) 2 + |A| 2 



J"B <H0 



y"B >(*{)■ 
(99) 

Thus the pions and the sigma meson will undergo a Mott tran- 
sition when their masses become larger than the two-particle 
continuum (Oq q , i.e., // B > a> qq . Using the mean-field results 
for A and M, we can calculate the two-particle continuum co qq 
as a function of /u B , which is shown in Fig|9] We find that the 
Mott transition does occur at a chemical potential fi B = //mi 
which is sensitive to the value of M*. The values of //mi for 
the four model parameter sets are shown in TablellVI For rea- 
sonable model parameter sets, the value of /(mi is in the range 
(7 - I0)m„. Above this chemical potential, the mesons are no 
longer stable bound states and can decay into quark-antiquark 
pairs even at zero momentum. We note that the Mott transition 
takes place well above the chiral restoration, in contrast to the 
pure finite temperature case where the mesons are dissociated 
once the chiral symmetry is restored 1841 [8511 . 

On the other hand, we find from the explicit forms of the 
meson propagators in AppendixlBlthat the decay process ji — ♦ 
qq is also possible at q + (even though |q| is small) due to the 
presence of superfluidity. Thus, we have another unusual Mott 
transition in the superfluid phase. Notice that this process is 
not in contradiction to the baryon number conservation law, 
since the Ub(1) baryon number symmetry is spontaneously 
broken in the superfluid phase. Quantitatively, this transition 
occurs when the meson mass becomes larger than the two- 
particle continuum co qq for the decay process n — > qq at q = 
+ . In this case, we have 



0)r. 



= J2 A /(M-f) 2 + |A| 2 , // B <// 
( 2|A| , // B > no- 



(100) 



The two-particle continuum a> qq is also shown in Figj9] We 
find that the unusual Mott transition does occur at another 
chemical potential // B = //M2 which is also sensitive to the 
value of M*. The values of /(M2 for the four model param- 
eter sets are also shown in TablellVI For reasonable model 
parameter sets, this value is in the range (5 - 8)m„. This pro- 
cess can also occur in the 2SC phase of quark matter in the 
N c — 3 case 118711 . In the 2SC phase, the symmetry break- 
ing pattern is SU C (3)®U B (1) ->SU C (2) ® U B (1) where the 
generator of the residue baryon number symmetry U B (1) is 
B = B — 2rg/V3 = diag(0, 0, 1) corresponding to the un- 



Set 


i 


2 


3 


4 


//mi [m K ] 


7.22 


7.76 


8.63 


9.62 


mm2 [»y 


5.29 


6.06 


6.96 


7.92 



TABLE IV: The chemical potentials yU M1 and ji M2 (in units of m n ) for 
different model parameter sets. 



paired blue quarks. Thus the baryon number symmetry for 
the paired red and green quarks are broken and our results can 
be applied. To show this explicitly, we write down the explicit 
form of the polarization function for pions in the 2SC phase 



nf c (g) 



n 2-coi ol(0 + 2 Tr[& (K)i 75 g (P)iy 5 ], (101) 



where Qa{K) is the propagator for the unpaired blue quarks. 
Here nl~ col ° r (Q) is given by (l29b (the effective quarks mass M 
and the pairing gap A should be given by the N c = 3 case of 
course) and corresponds to the contribution from the paired 
red and green sectors. The second term is the contribution 
from the unpaired blue quarks. Therefore, the unusual decay 
process is only available for the paired quarks. 



V. BEYOND-MEAN-FIELD CORRECTIONS 

The investigations in Sec. ITTT1 and llVl are restricted in the 
mean-field approximation, even though the bosonic collec- 
tive excitations are studied. In this section, we will include 
the Gaussian fluctuations in the thermodynamic potential, and 
thus really go beyond the mean field. The scheme of go- 
ing beyond the mean field is somewhat like those done in 
the study of finite temperature thermodynamics of the NJL 
model 1881 18911 ; however, in this paper we will focus on the 
beyond-mean-field corrections at zero temperature, i.e., the 
pure quantum fluctuations. We will first derive the thermo- 
dynamic potential beyond the mean field which is valid at ar- 
bitrary chemical potential and temperature, and then briefly 
discuss the beyond-mean-field corrections near the quantum 
phase transition. The numerical calculations are deferred for 
future studies. 



A. Thermodynamic potential beyond the mean field 

In the Gaussian approximation, the partition function can 
be expressed as 

Z N jl - exp(-<S<° ff ) ) f[do-][dn][dtf][d(/>] exp( - 5^).(102) 

Integrating out the Gaussian fluctuations, we can express the 
total thermodynamic potential as 



n(T,n B ) = n (T, f i B ) + si a (T,ij B ), 



(103) 



where the contribution from the Gaussian fluctuations can be 
written as 

Q fl = - Y[lndetM(0 + lndetN(0]. (104) 



However, there is a problem with the above expression, 
since it is actually ill-defined: the sum over the boson Mat- 
subara frequency is divergent and we need appropriate con- 
vergent factors to make it meaningful. In the simpler case 
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without superfluidity, the convergent factor is simply given by 
e 'v„fi* |gg[[g9}i i n our case, the situation is somewhat different 
due to the introduction of the Nambu-Gor'kov spinors. Keep 
in mind that in the equal time limit, there are additional factors 
e'"» 0+ for @u(K) and <r''<"» 0+ for @ 12 {K). Therefore, to get the 
proper convergent factors for Qfl, we should keep these factors 
when we make the sum over the fermion Matsubara frequency 
cj n in evaluating the polarization functions rLj(0 and 11^(0. 
The problem in the expression of Qfl is thus from the op- 
posite convergent factors for My and M22- From the above 
arguments, we find that there is a factor e' v, "° for Mu and 
e -'v„,o f or jy[ 22 Keep in mind that the Matsubara sum £ m is 
converted to a standard contour integral (iv m — > z). The con- 
vergence for z — > +00 is automatically guaranteed by the Bose 
distribution function b{z) = \j(eP z - 1), we thus should treat 
only the problem for z — > — °°. To this end, we write the first 
term of Qfl as 



^lndetM(0 

e 



Xi[ lnM i 



lnM22<? 



-iv,„0 + 



+ ln 



detM 



,MllM 22/ 

Using the fact that M 22 (0 = Mn(-g), we obtain 



^lndetM(0 



Z 1 " 

Q 



Mu(fi) 



M 22 (0 



detM(0 



(105) 
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Therefore, the well-defined form of Qfl is given by the above 
formula together with the other term 2glndetN(0 associ- 
ated with a factor e n '"'° . 

The Matsubara sum can be written as the contour integral 
via the theorem 2mgO' v »i) = § c dz/(2m)b(z)g(z), where C 
runs on either side of the imaginary z axis, enclosing it coun- 
terclockwise. Distorting the contour to run above and below 
the real axis, we obtain 

-<5 22 (w,q) + 3^(w,q)], (107) 

where the scattering phases are defined as 

<!>m(w, q) = ImlndetM(w + /0 + ,q), 

6u(oJ,q) = ImlnMn(w + /0 + ,q), 

5ii((x), q) = Im In M22(w + /0 + , q), 

S„(cj, q) = Imln [(2G) _1 + U n (cu + /0 + , q)] . (108) 

Keep in mind the pressure of the vacuum should be zero, the 
physical thermodynamic potential at finite temperature and 
chemical potential should be defined as 

Q p hy(T ) /iB) = n(r, J u B )- 0(0,0). (109) 



B. Thermodynamic consistency of the vacuum 

As we have shown in the mean-field theory, at T = 0, the 
vacuum state is restricted in the region |yU B | < m n . In this re- 



gion, all thermodynamic quantities should keep zero, no mat- 
ter how large the value of fig is. While this should be an obvi- 
ous physical conclusion, it is important to check whether our 
beyond-mean-field theory satisfies this condition. 

Notice that the physical thermodynamic potential is defined 
as Qphy(^B) = Q(^b) _ Q(0), we therefore should prove that 
the thermodynamic potential Q(^b) keeps a constant in the 
region \jib\ < m„. For the mean-field part Qo, the proof is 
quite easy. Because of the fact that M* > m„/2, the solution 
for M is always given by M = M*. Thus Qo keeps its value at 
fiB - in the region |//bI < m n . 

Now we turn to the complicated part Qfl. Since A = 0, all 
the off-diagonal elements of M vanishes, and Qfl is reduced to 



Q1-1 



15> 



h + ^ {Q) 



2G +U * (Q) 



Jv m 0* 



Jv,„a* 



J> 



— + n d (0 

4G 



(110) 



where Yi a {Q) - 1133(0 and we should set A = and M - M„ 
in evaluating the polarization functions. First, we can easily 
show that the contributions from the sigma meson and pions 
do not have explicit yU B dependence and thus keep the same 
values as those at /i B = 0. In fact, since the effective quark 
mass M keeps its vacuum value M, guaranteed by the mean- 
field part, all the fis dependence in 11^^(0 is included in the 
Fermi distribution functions f(E ± fis/2). Since M, > /j.b/2, 
they vanish automatically at T - 0. In fact, from the explicit 
expressions for 11,^(0 in Appendix [B] we can check that 
there is noyUB independence in 11^^(0. 

The diquark contribution, however, has an explicit yU B de- 
pendence through the combination iv m +/u B in the polarization 
function IId(0. The diquark contribution (at T — 0) can be 
written as 

vn f ° dOJ 

Qd = - V — S d (a>, q), 

5 d (w,q) = Imln[(4Gr I +n d (w + /0 + ,q)]. (Ill) 

Making a shift co — > ai - /^b, and noticing that fact U.d(a> - 
jUb, q) = n^(w, q)/2, we obtain 



Qh = -Y 



q 



— tf„(w,q). (112) 

J-co 7T 



To show the above quantity is in fact /i B independent, we sep- 
arate it into a pole part and a continuum part. There is a 
well-defined two-particle continuum E c (q) for pions at arbi- 
trary momentum q, 



E c (q) = min k (#* + £* +q ) . 



(113) 



The pion propagator has two symmetric poles ±a> n (q) when q 
satisfies «„-(q) < E c (q). Thus in the region \a>\ < E c (q), the 
scattering phase 8„ can be analytically evaluated as 

8 n (u), q) = n [0 (-oj - aj n (q)) -®(oj- w„(q))] . (1 14) 
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Since E c (q) > b)„(q) > m„ > fis, the thermodynamic potential 
Qd can be separated as 

" d = 2 K(q) - E M)] - X f " ~ 6n(0J ' qX (1 15) 

which is indeed /j.b independent. Notice that in the first term 
the integral over q is restricted in the region |q| < q c where q c 
is defined as oj^iqc) — E c (q c ). 

In conclusion, we have shown that the thermodynamic po- 
tential Q. in the Gaussian approximation keeps a constant in 
the vacuum state, i.e., at |/jbI < m n and at T — 0. All other 
thermodynamic quantities such as the baryon number density 
keep zero in the vacuum. The subtraction term Q(0, 0) in the 
Gaussian approximation can be expressed as 

G(0, 0) = Q vac (M,) + | £ K(q) - E c (q)] 



X-£.,-iqj j 
— [tfo.(w,q) + 5^(ca,q)]. (1 



16) 



C. Quantum corrections near the phase transition 

Now we consider the beyond-mean-field corrections near 
the quantum phase transition point hb — m n . Notice that the 
effective quark mass M and the diquark condensate A are de- 
termined at the mean-field level, and the beyond-mean-field 
corrections are possible only through the equations of state. 

Formally, the Gaussian contribution to the thermodynamic 
potential Qfl is a function of hb,M and y - |A| 2 , i.e., Qfl = 
Qn(/JB,y, M). In the superfluid phase, the total baryon density 
including the Gaussian contribution can be evaluated as 



n(jj B ) = Hoi/An) + waO"b), 



(117) 



where the mean-field part is simply given by «o(//b) = 
-dtlo/d/iB and the Gaussian contribution can be expressed as 



«h(/Ub) = - 



dQ fl <9Q fl dy 3Q fl dM 
dfiB dy d/jB dM dfiB 



(118) 



The physical values of M and |A| 2 should be determined by 
their mean-field gap equations. In fact, from the gap equations 
dtlo/dM = and dQ^/dy - 0, we obtain 



<9 2 Q <3 2 Q dy d 2 Q dM 



+ 



+ 



dfindM ' dydM dfi B dM 2 d/i B 
<9 2 Q <3 2 Q) & d2Q o dM 



0. 



dfiftdy dy 2 dfiB dMdy dfis 



= 0. 



(119) 



Thus, we can obtain the derivatives dM/dfiB and dy/dfiB ana- 
lytically. Finally, «fl(/iB) is a continuous function of /j.b guar- 
anteed by the properties of second order phase transition, and 
we have «a(m ff ) = 0. 

Next we focus on the beyond-mean-field corrections near 
the quantum phase transition. Since the diquark condensate 
A is vanishingly small, we can expand the Gaussian part Qfl 



in powers of |A| 2 . Notice that jjb and M can be evaluated 
as functions of |A| 2 from the Ginzburg-Landau potential and 
mean-field gap equations. Thus to order (9(|A| 2 ), the expan- 
sion takes the form 



Q fl - ^|A| 2 , 

where the expansion coefficient 77 is defined as 



n = 



dClft dQfl du.B <9Qfl dM 

-1 •— + 

dy dfiB dy dM dy 



(120) 



(121) 



fi B =m„y=0,M= Af, 



Using the definition of rag, we find that 77 can be related to «g 

by 



J] - mtmn)— — 
ay 



y=0 



(122) 



Thus, the coefficient 77 vanishes, and the leading order of the 
expansion should be <9(|A| 4 ). 

As shown above, to leading order, the expansion of Qfl can 
be formally expressed as 



Qfl 



-fySIAI 4 . 



(123) 



The method to derive the exact expression of the numerical 
factor ( is shown in Appendix|C] Notice that the factor ( is in 
fact hb independent, thus the total baryon density to leading 
order is 



n = n + fiS|A| 



2 d|A| 2 



d/iB 



(124) 



Mb="% 



Near the quantum phase transition point, the mean-field con- 
tribution is «o = |\|*ol 2 = 2m„-J"|A| 2 from the Gross-Pitaevskii 
free energy. The last term can be evaluated using the analyti- 
cal result 



1 a v& 
Ivor = — 

go 



|A|" 



2m„S 



-f*d, 



(125) 



which is in fact the solution of the mean-field gap equations. 
Therefore, to leading order, the total baryon density reads 

n=(l+02m n J\A\ 2 . (126) 

On the other hand, the total pressure P can be expressed as 



P = (1+^|A| 4 . 



(127) 



Thus we find that the leading order quantum corrections are 
totally included in the numerical factor £ . Setting ( — 0, we 
recover the mean-field results obtained in Sec. [Hi] 

Including the quantum fluctuations, the equations of state 
shown in d65ll are modified to be 



P(n) = 
fi 3 (n) = m 



1 2na dd 



1 + ( m„ 

1 47ra dd 

; 

1 + C m„ 



(128) 
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This means, to leading order, the effect of quantum fluctua- 
tions is giving a correction to the diquark-diquark scattering 
length. The renormalized scattering length is 



J dd 



«dd 

T+7' 



(129) 



Generally, we have ( > and the renormalized scattering 
length is smaller than the mean-field result. 

An exact calculation of the numerical factor £ can be per- 
formed using the method shown in Appendix[C] However, this 
needs huge numerical power and we defer it to future work 
119011 . In this paper we will give an analytical estimation of £ 
based on the fact that the quantum fluctuations are dominated 
by the gapless Goldstone mode. To this end, we approximate 
the Gaussian contribution Qg as 

Qfl * \j] ln [ D ^ (Q)D ^ ( ' Q) + 3 ^ A]4 



+2y3|A| 2 (D d \Q) + D d l (-Q)) 



(130) 



where D d l (Q) is given by (l54l and can be approximated by 
( 1561 . Subtracting the value of Qfl at /i B = m n with A = 0, and 
using the result jU B — m^ + golYol 2 from the Gross-Pitaevskii 
equation, we find that (, can be evaluated as 



B m 

<r = ^(/i + / 2 )--f(/i + /2), 

<J Jn 

where the numerical factors l\ and h are given by 



Z 2 + X 2 



(131) 



2 Zj Zj (Z 2 - X 2 ) 2 - 4Z 2 

m X '" m 



h 



4 ZZ 

m X i(Z m 



(3Z 2 - X 2 ) 



2^2 



X 2 ) 2 -4Z 2 ,r 



(132) 



Here the dimensionless notations Z m and X are defined as 
Z m = iv m /m K and X = q/m n respectively. Notice that the 
integral over X is divergent and hence such an estimation has 
no prediction power due to the fact that the NJL model is non- 
renormalizable. However, regardless of the numerical factor 
I\ + I2, we find that f oc ra 2 //, 2 . Thus, the correction should 
be small in the nonlinear sigma model limit m„ «; 2M*. 



The transition temperature T c is determined by the Thouless 
criterion D d '(0, 0) = which can be shown to be consistent 
with the saddle point condition 6S e ff/6(f>\^=o = 0. Its explicit 
form is a BCS-type gap equation 



45=^'?/ 



d\ l-2/(ff) 
(2tt)3 2£ ■ 



(133) 



Meanwhile, the dynamic quark mass M satisfies the mean- 
field gap equation 



M - mo 
2GM 



r d'k 



k !-/(£)-/(£) 



£ k 



(134) 



To obtain the transition temperature as a function of 
n, we need the so-called number equation given by n - 
-dQ./d/j.B, which includes both the mean-field contribution 

«o(j«b, T) = 2Nj 2k /(£k) _ /(^k) anc ' tne G auss i an contri- 
bution «fl(yUB, T) = -dQ.fi/diJB- At the transition temperature 
where A = 0, Qfl can be expressed as 

x [26 d (co,q) + Ww,q) + 36 n (o,q)], (135) 

where the scattering phases are defined as 6a{<jd, q) = 
Imln[l/(4G) + U d (aj + z'0 + ,q)] for the diquarks, 6a-(a>,q) = 
Imln[l/(2G) + ^(w + ;0 + ,q)] for the sigma meson and 
6„(u),q) = Imln[l/(2G) + Ii n {u) + /0 + ,q)] for the pions. Ob- 
viously, the polarization functions should take their forms at 
finite temperature in the normal phase. 

The transition temperature T c at arbitrary baryon number 
density n can be determined numerically via solving simulta- 
neously the gap and number equations. However, in the dilute 
limit « — > which we are interested in this section, analyti- 
cal result can be achieved. Keep in mind that T c — > when 
n — > 0, we find that the Fermi distribution functions /(^) 
vanish exponentially (since M, - m n j2 » T c ) and we ob- 
tain yUB = m„ and M — M* from the gap Eqs. d 133b and 
dl34l ), respectively. Meanwhile the mean-field contribution of 
the density «o can be neglected and the total density n is thus 
dominated by the Gaussian part «fl. When T — > we can 
show that Yla-ico, q) and Yl n {co, q) are independent of fig, and 
the number equation is reduced to 



D. Transition temperature 

While the effect of the Gaussian fluctuations at zero temper- 
ature is to give a small correction to the diquark-diquark scat- 
tering length and the equations of state, it can be significant at 
finite temperature. In fact, as the temperature approaches the 
critical value of superfluidity, the Gaussian fluctuations should 
dominate. In this part, we will show that to get a correct crit- 
ical temperature in terms of the baryon density n, we must 
go beyond the mean field. The situation is analogous to the 
Nozieres-Schmitt-Rink treatment of molecular condensation 
in strongly interacting Fermi gases 



Zf °° da 



doj <95 d (w,q) 
b(aji)- 



<3yU B 



(136) 



Since T c — > 0, the inverse diquark propagator can be reduced 
to !D d l ((L>, q) in (l54l i. Thus the scattering phase 6d can be well 
approximated by 5 d (u>, q) = 7r[0(// B -6 q -w)-0(w-;U B -eq)] 
with e q = sjq 2 + m 2 . Therefore, the number equation can be 
further reduced to the well-known equation for ideal Bose- 
Einstein condensation, 



2j [ b ( € <i - J«b) - b(e q + Mb)] 



(137) 
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Since the above equation is valid only in the low density limit 
n — * 0, the critical temperature is thus given by the nonrela- 
tivistic result 



2n 

T r = — 



£(3/2) 



2/3 



(138) 



At finite density but «a^ d <k 1, there exists a correction to T c 



which is proportional to n^ 3 addll3i|. Such a correction is hard 
to handle analytically in our model since we should consider 
simultaneously the corrections to M and fig, as well as the 
contribution from the sigma meson and pions. 



VI. SUMMARY 

In summary, we have examined the NJL model descrip- 
tion of weakly interacting Bose condensate and BEC-BCS 
crossover in QCD-like theories at finite baryon density. Our 
main conclusions are as follows: 

(l)Near the quantum phase transition point /iB = m n , we have 
performed a Ginzburg-Landau expansion of the effective po- 
tential. At the mean-field level, the Ginzburg-Landau free en- 
ergy is essentially the Gross-Pitaevskii free energy describ- 
ing weakly repulsive Bose condensates after a proper redefi- 
nition of the condensate wave function. The obtained diquark- 
diquark scattering length reads Add = tn n l{\6nf%), which re- 
covers the tree-level result predicted by chiral Lagrangian. 
(2)We have analytically shown that the Goldstone mode near 
the quantum phase transition point takes the same disper- 
sion as the Bogoliubov excitation in weakly interacting Bose 
condensates, which gives a diquark-diquark scattering length 
identical to that in the Gross-Pitaevskii free energy. The mix- 
ing between the sigma meson and the diquarks plays an im- 
portant role in recovering the Bogoliubov dispersion. 
(3)The results of baryon number density and in-medium chi- 
ral and diquark condensates predicted by chiral perturbation 
theory are analytically recovered near the quantum phase tran- 
sition point in the NJL model. 

(4)At high density, the superfluid matter undergoes a BEC- 
BCS crossover at /jb - (tn a /m n ) l ^m„ =* (1.6 - 2)m JT . At 
I-Lb - 3ntn, the chiral symmetry is approximated restored and 
the spectra of pions and sigma meson become nearly degen- 
erate. Well above the chemical potential of chiral symmetry 
restoration, the degenerate pions and sigma meson undergo a 
Mott transition, where they become unstable resonances. Be- 
cause of the spontaneous breaking of baryon number sym- 
metry, mesons can decay into quark pairs in the superfluid 
medium at nonzero momentum. 

(5)The general theoretical framework of the thermodynamics 
beyond the mean field is established. It is shown that the vac- 
uum state in the region |//bI < m„ is thermodynamically con- 
sistent in the Gaussian approximation, i.e., all thermodynamic 
quantities keep vanishing for |//bI < tn„ even though the Gaus- 
sian fluctuations are included. 

(6)Near the quantum phase transition point, we find that the 
effect of the leading order beyond-mean-field correction is to 
renormalize the diquark-diquark scattering length. The cor- 
rection to the mean-field result is estimated to be proportional 



to m 2 n l f-. Our theoretical approach provides a new way to cal- 
culate the diquark-diquark or meson-meson scattering lengths 
in the NJL model beyond the mean-field approximation. We 
also find that we can obtain a correct transition temperature of 
Bose condensation in the dilute limit once the beyond-mean- 
field corrections are included. 

Our studies can be generalized to describepion conden- 
sation at finite isospin chemical potential [i\ |4], |5[] and kaon 
condensation at finite strangeness chemical potential /is 19111 ■ 
In the NJL model, pion condensation is shown to occur at 
|/ii| = m„ when |//s| < 'Wk - m„/2, and kaon is shown to 
condense at \fj.s\ - m^ when ^i\ — > J46l 14811 . The general- 
ization to pion condensation is straightforward. The obtained 
Ginzburg-Landau and Gross-Pitaevskii free energies are the 
same as those derived in this paper, if we replace //b —* Hi- 
The results are valid both for N c = 2 and N c = 3 cases. At 
the mean-field level, the results for diquark condensation at 
N c — 2 and pion condensation at N c — 3 are formally iden- 
tical in the NJL model. Significant difference may appear 
if we consider the beyond-mean-field corrections, since for 
the N c = 3 case the scalar diquarks are not pseudo-Goldstone 
bosons. A calculation of the pion-pion scattering length in the 
1-2 channel can be performed within our theoretical frame- 
work. The calculations of kaon condensation and kaon-kaon 
scattering length are also possible, but somewhat complicated 
due to the large mass difference between the light and strange 
quarks. 

We are also interested in how the beyond-mean-field cor- 
rections modify the superfluid equations of state. As we learn 
from the knowledge of BEC-BCS crossover in cold Fermi 
gases, the superfluid equations of state can be strongly mod- 
ified in the crossover regime 11571 15811 . corresponding to the 
moderate baryon density in our case. This issue is also im- 
portant to the color-superconducting quark matter I92U94I1 at 
moderate density, i.e., for quark chemical potential around 
400MeV where the pairing gap can be of order (9(100MeV). 
The numerical works are in progress. 
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Appendix A: Fermionic Model Description of Dilute Bose 
Condensate 



In this appendix, we briefly review the theory of molecu- 
lar Bose condensation in two-component Fermi gases in the 
strong coupling limit. While there exist many theoretical ap- 
proaches J95I497I1 to deal with this problem, we employ the 
field theoretical approach 11571 15811 parallel to that used in this 
paper. 

The Lagrangian density of the system can be written as 

L= E^r + 2^ + T <7+g ^^ r (A1) 
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where i//^ denote the two-component (nonrelativistic) 
fermion fields with equal masses m and chemical potentials [i. 
The gas is assumed to be dilute, and the coupling constant g 
can be related to the s-wave fermion-fermion scattering length 
a s as 



Ana. 



— V — 

's( A ) + ^2e k ' 



(A2) 



where e k = k 2 /(2m). In the dilute limit, we can take the limit 
A — > oo in the final result. 

Performing the Hubbard-Stratonovich transformation with 
the auxiliary boson field <p(x) - gifr.(x)iffJx), and defin- 
ing the Nambu-Gor'kov representation T*' = (iff., iff.), we 
can evaluate the partition function of the system as Z = 
$[dW>][cT¥][d<p>][d(f>] expC-^ff), where 



lA 



elf 



fdx-^— -- fdx f dx /y frt(x)G- l (x,x'y¥(xXA3) 



and the inverse fermion propagator G ' is given by 



4>\x) -d T 



2iii 



5(x - x). 



(A4) 



Then integrating out the fermionic degree of freedom, we get 
Z = f[d(fi][d<f>] exp(-»Seff) where the bosonized effective ac- 
tion reads 



•Seff [</> , 4>] = dx 



/• 



\<P(x)\ 2 



TrlnG" 1 ^,*')- 



(A5) 



1. Mean-field theory 

First, we consider the mean-field theory where the auxiliary 
boson field <f>(x) is replaced by its expectation value {<p(x)) — 
A(x). In the strong coupling limit a s — > + , the fermion chem- 
ical potential fi approaches --E/,/2 with E\, - 1 /(ma 2 ) being 
the molecular binding energy. Since the pairing gap |A| «: \fi\, 
we can expand the effective action in powers of |A|, which 
resulting in a Ginzburg-Landau free energy functional 

Vgl[A(x)] = fdx aW*^ -7V 2 )a(x) 



+ a|A(x)| 2 + -/?|A(x)f 



(A6) 



The coefficients a, ft of the potential terms can be ob- 
tained from the mean-field thermodynamic potential Qo = 
(T/V)S e ff[A i , A] which can be evaluated as 



Oo 



Ana. 



|A| 2 



ltf- S *- 6 -£ 



(A7) 



where £k = «k - H and E^ - Jj£ + |A| 2 . After a simple 



algebra, the coefficients a and /? can be evaluated as 



1 



= J- n I V"2"V ' 



P 



1 



8tt (-Imitfl 1 ' 



(A8) 



From the expression of a, we find that a quantum phase 
transition from vacuum to Bose condensation takes place at 
li - -l/(2ma 2 ) = —Eb/2. Thus, near the phase transition, a 
can be simplified as 



m 2 a s 



a =a -- 



8/r 



■fJb, 



(A9) 



where ii b — 2ii + Ej, is the boson chemical potential. Further, 
setting /u - —Eb/2, /3 can be simplified as 



P 



\6n 



(A10) 



The coefficients y, k of the kinetic terms can be obtained 
from the inverse boson propagator !D~ l (Q) with A = 0. It can 
be evaluated as 



O l (Q) = - 



Ana s 



z 



1 



1 



■ £k - £k+q 2e k 



(All) 



In the strong coupling limit, it can be well approximated 
asS 



£T'(0-- 



8/r 



IVm- 



Am 



(A12) 



In summary, if we define the new condensate wave function 
y(x) by 



y(x) = 



8/r 



AW, 



(A13) 



the Ginzburg-Landau free energy can be reduced to the Gross- 
Pitaevskii free energy of dilute Bose gases, 



VgpLVW] 



fdxy\x)i 



d V 2 \ 

v,(x)| ^-2^ v(x) 



2 1 Ana bh 4 

tx b \v(x)r + - kWl 

2 m b 



(A14) 



where m b - 2m is the boson mass and a bb = 2a s is the 
boson-boson scattering length. Since a s — > + , the interac- 
tions among the composite bosons are repulsive and weak. 



2. Beyond-mean-field corrections 

To study the beyond-mean-field corrections, we consider 
the fluctuations around the mean field. Making the field shift 
<f>(x) — > A + <p(x), we can expand the effective action S e s in 
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powers of the fluctuations. The zeroth order term S~2 is just 
the mean-field result, and the linear terms vanish automati- 
cally guaranteed by the saddle point condition for A. The 
quadratic terms, corresponding to Gaussian fluctuations, can 
be evaluated as 

^ = ^(^)«-e))M(2)(;S), (M5) 

where the inverse boson propagator M is given by 

Mn (g) = MzsC-fi) = - + Y622(K)g n (K + Q) 



-i-Z 



«k"Lq 



2 2 
k k+q 



/V,„ - £k - £k+q 'V„, + £k + £k 



+q 



(A16) 



and 



Mi 2 (0 = M 2 i(<2) = J] Gii(K)@ n (K + 



MkVkMk+q^k+q 



Z; Kk'kKk+q 1 
\ iv m + E k + 



MkVk«k+qVk+q 



J"* \ iv m + £k + ^k+q 'V,,, - E\s. - £k+q 



(A17) 



Here the fermion Green function Q is defined as Q ' = 
G ' [A] and the BCS distribution functions v 2 = (1 -£ k /E k )/2 
and m^ = 1 - v 2 . are used. 

In the strong coupling limit where |A|/|/i| <*: 1, the matrix 
elements of M can be analytically evaluated. We have 115711 



M(0 



m 2 a, 



8/r 



~ M + 2^ " A*b + 2golVol 2 

golvol 2 



golVol 2 

W + 2^ _ ^ + 2 ^ol¥o| : 



(A18) 



where go = Ana^/m^ and |yo| 2 is the minimum of the Gross- 
Pitaevskii free energy. Together with the mean-field result for 
the boson density ri\, = |yol 2 , we can show that the Goldstone 
mode takes a dispersion relation given by 



w(q) 



/ q 2 / q 2 87ra b b»b 
I2mb \2mb mb 



(A19) 



which is just the Bogoliubov excitation in a dilute Bose con- 
densate. 

To evaluate the thermodynamic potential beyond the mean 
field, we express the partition function in the Gaussian ap- 
proximation as 

Z - exp (-SQ JW] [#] exp ( - S ( ^). (A20) 

Integrating out the Gaussian fluctuations, the total thermody- 
namic potential can be expressed as 



QOu) = fioQz) + QflOu), 



(A21) 



where the contribution from the Gaussian fluctuations can be 
evaluated as 115811 



Q fl 



&* 



Mn(fi) 



M 22 (0 



detM(0 



(A22) 



Near the quantum phase transition point fi = -Et,/2, we 
can expand Qh in powers of |A| 2 . Because of the properties 
of second order phase transition, the terms of order <9(|A| 2 ) 
vanish. To leading order, the result is 15811 
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256/r 
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where the numerical factor £ - 2.61. From the Gross- 
Pitaevskii free energy, we find that the pressure P in the mean- 
field approximation can be expressed as 



3 3 



32tt 



W 



(A24) 



Thus, to leading order, the beyond-mean-field corrections 
renormalize the boson-boson scattering length Abb- The new 
renormalized scattering length reads 



«bb 



2a s 

T+7 



0.55a 5 



(A25) 



Notice that this result is quite close to the exact result for the 
four body problem of 0.6a s 19811 . This means the quantum 
fluctuations are almost correctly included in the present theo- 
retical approach. 

Further, going beyond the leading order we find that we can 
fit Qfl to the functional form 115811 



" fl = 2a? ( Cl ^ + C2 ^ /2 + " ' ^ ' 



(A26) 



where /5b = Hb/Eb and the dimensionless factors c\, c 2 can be 
numerically determined. Solving for the molecular chemical 
potential //b one obtains 



^b = 



Anavbnb 
m b 



32 / 3 xi 

1 + f^-p ("bflbbj 

3 V 7r 



(A27) 



with the coefficient £, - 0.94|58[| which is 6% smaller than the 
Lee-Huang- Yang result £ = 1 IU5I1 . 



Appendix B: The One-Loop Susceptibilities 

In this appendix, we evaluate the explicit forms of the one- 
loop susceptibilities n;j(0 (i,j = 1,2,3) and n^(0. At arbi- 
trary temperature, their expressions are rather huge. However, 
at T = 0, they can be written in rather compact forms. For 
convenience, we define A = \A\e' e in this appendix. 

(I) Diquark sector. First, the polarization functions Tl n (Q) 
and ni 2 (0 can be evaluated as 



26 



IIii(0 = KN^ 



k L 



k L 

+ 



(« k ) 2 (« B ) 2 



(« k ) 2 (vp~) 2 

'm — ^k — P 



(v k ) 2 (v D ) 2 



«) 2 K) 2 



niaCfi) = N c NfY, ( 



i iv m + E k + E v 
"k v k"p v p 



/v m + E k + E p 


(v k ) 2 («^) 2 


iv„, + E k + £+ 


«4 v k«j v ; 


^Vm — ^^ — &p 


\\K v p 



iv,„ + £++£+ 


«) 2 (vp 2 


iv,„ + £+ + £ p 


M k v k M P v P 


i'v m +£+ + £■+ 


M k V+ Mp Vp 



(v^) 2 (vp 2 

(<) 2 (« D ) 2 



r + 



r. 



iv,„ + E. 



iv,„ +E++E p 



Wm ^jj ^p 

M k v k M P y P 

iv m - E+- E p 



p ) 



T. 



,-•« 



(Bl) 
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where p = k + q. Here T± are factors arising from the trace in 
spin space, 

T ± = U k Z + f, (B2) 

2 2EkEp 

and ut, v k are the BCS distribution functions defined as 



Thus, near the quantum phase transition point, we have 
rii2(0) = A 2 /?] + <9(|A| 4 ). On the other hand, a simple algebra 
shows that 



4G +IIll(0) 



111,2(0)1 
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At Q - 0, we find that 

n 12 (0) = A 2 jN c N f J] 



(£ k ) 3 



(Et) 3 



(B3) Therefore, the mean-field gap equation for A ensures the 
Goldstone's theorem in the superfluid phase. 

(II) Diquark-sigma mixing terms. The term II13 standing 
for the mixing between the sigma meson and the diquarks 
(B4) reads 
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where the factors 1+ are defined as 



Ml 1 1 
J ± = - — ± — 



2 \£ k 



(B7) 



One can easily find that II13 ~ MA, thus it vanishes when A 
or M approaches zero. At Q — 0, we have 



1 v-i M 

n 13 (0) = A-iv c iv>V — 

^ 1, ^k 
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Thus the quantity Hq defined in ( TTTb can be evaluated as 

1 ... vnxnM, 1 



H a 



e=± k 
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(III) Sigma meson andpions. The polarization function II33 
which stands for the sigma meson can be evaluated as 
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where the factors T' ± are defined as 

lk_p^ (BU) 

2 2£ k £ , p 

At 2 = and for A = 0, we find that 

1 n 1 vnM? 

M 33 (0) = — - 2^ £ - + 2iV c ^ J] ^ 
5 2 n (y, M) 



k £ * 



k k 



<9M 2 



(B12) 



y=0 



Finally, the polarization function Yi n (Q) for pions can be ob- 
tained by replacing T' ± — > 7"±. Thus, when M — > 0, the sigma 
meson and pions become degenerate and chiral symmetry is 
restored. 



Appendix C: Expansion of Q. a in Terms of |A| 2 

In this appendix, we derive the expression of the Taylor ex- 
pansion of Qfl in terms of |A| 2 = y. As we have shown in Sec. 
|VJ the leading-order term should be <9(|A| 4 ). Thus, we need 
to evaluate the numerical factor f , A key problem here is that 
the effective quark mass M and the chemical potential fis are 
both functions of |A| 2 determined at the mean-field level. 

First, we expand the matrix elements of M and N in terms 
of y. Any of these elements denoted by F is a function of 



J 



Us, M and y. Our method of expansion is as follows. We 
firstly expand F(jjb , M, y) in terms of y formally with //b and 
M being fixed parameters, i.e., 

F(jiB,M,y) = Fo(fiB,M) + F l (fiB,M)y + F 2 (HB,M)y 2 + 0(y 3 ), 

(CI) 
where Fo(/ib,M) = F(jib, M, 0) and the expansion coefficients 
are defined as 



F 1 ( f i B ,M) = 



dFQto,M,y) 



dy 



F 2 (MB,M) = - — 2 



v=0 



(C2) 



We then expand the coefficients Fi(jib,M) (i = 0,1,2) at 
(jU B , M) - (m n , M t ), using the fact that 



1 
M^M t v, 



Hb - m„ + 



Doing this we formally obtain 



2m n J 



y- 



- i 
Finally, up to order 0(y 2 ), we have 



+ F 2 (m n ,M*)y 2 + 0(y 3 ). 



(C3) 



(C4) 



F(jx B , M,y) = F (m„, M.) + [F l Q (m n , M t ) + F y (m„, M,)\y + [pfon*, M.) + F\(m n , M.) + F 2 (m„, M.)] 

Using this method, we can expand the matrix elements of M and N formally as follows: 

Mn(S) = D' d -\Q) + Xi(fi)|A| 2 + T,(0|A| 4 + 0(|A| 6 ), 
M 12 (0 = A 2 Z(0 + 0(|A| 4 ), 
M 13 (0 = AW(Q) + 0(|A| 3 ), 

M 33 (0 = D* o r\Q)+X 2 (Q)\A\ 2 + Y 2 (Q)\A\ 4 + 0(|A| 6 ), 
Nn(0 = D;-\Q) + X i (Q)\A\ 2 + F 3 (0|A| 4 + 0(|A| 6 ). 



v 2 . 



(C5) 



(C6) 



Here D* d l {Q) is defined as D d \Q;^ B = m n ). 
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Meanwhile, the thermodynamic potential Qfl can be expressed as 



Q fl4z{ ln 

Q { 
Using the expansion 
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detM(Q) 
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we find that the factor f is given by 
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It is obvious that Z(0 = B(g) and W(0 = //(G) where fl(g) 
and //(G) are defined in (TTTT ). On the other hand, since 



/S 



1 



2m„J" 2M t 2M« 



(C9) 



to leading order of 0(m„/2M*), we can set jjb = »!* in all 
equations and identify Xi(0 = A(Q) defined in ( fTTT i. 
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